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ABSTRACT 


The FSCBG computer program described in this report combines and 
implements mathemattcal models for atreraft wake effects, line-source dis- 
perston, drop evaporation and canopy penetration. The computer program ts 
destgned to account for the atmospherte disperston, transport and depost- 
tton of all aertal spray material from the time of release unttl all spray 
matertal ts etther depostted or, in the case of spray drift beyond a target 
area, spray concentration and deposttion levels are insignificant. Spectfic 
calculations made by the computer program tnelude spray concentrations and 
dosages above forest canopies as well as the spray deposttion within and 
below forest canoptes resulting from aertal spray releases made along 
stngle or multtple flight paths. Important applications of the computer 
progran are the optimization of spray program destgn and operation wtth 
respect to the selection of atreraft spray systems, fltght altitudes, 
swath widths, spray rates, and scheduling of spray operations; evaluation 
and analysts of fteld measurements of spray deposition; and, assessments 


of the environmental impact or hazard posed by aertal spray operattons. 


Prepared for USDA-Forest Service, Forest Insect and Disease Management 
Application Group, Davis, CA under Contract No. 53-91S8-9-6127; John W. 
Barry, USDA-Forest Service, Project Officer. 
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SECTION 1 
INTRODUCTION 


Lee BACKGROUND 


The USDA Forest Service (USFS) uses aerial spray applications to 
prevent and suppress tree damage due to pest infestations. As part of an 
extensive program to develop and evaluate pest management systems, the 
USFS is interested in achieving a complete understanding of the behavior 
of spray material from the time spray is released from the aircraft until 
it has been deposited or, in the case of spray drift, dissipated to concen- 
tration/dosage levels that are environmentally insignificant. Because 
mathematical spray dispersion models are useful in determining interactions 
of the many factors affecting spray peer acioncn the USFS has supported 
the application and development of these models over the last decade. 
Simplified aerial line source models developed for the US Army were 
applied (GCA Corporation, 1971; H. E. Cramer Company, 1973) early in the 
decade to determine optimum swath widths and application rates for use in 
pilot tests of insecticides under consideration at that time for control 
applications in western forests. The implications of these early efforts 
in the use of mathematical models to improve the planning, conduct and 
analysis of spray programs were reviewed in a paper presented at the USFS 
Workshop for Aerial Application of Insecticides Against Forest Defoliators, 
held in Missoula, Montana, 23-29 April 1974 (Dumbauld, Cramer and Barry, 
1975). The H. E. Cramer Company (Dumbauld and Bjorklund, 1977) used aerial 
spray dispersion models to assist the State of Maine Bureau of Forestry 
in determining offset distances required for various aircraft to ensure 
that drift from spray blocks in the Maine 1977 spray program posed no 
environmental hazard to exclusion areas (waterways, homes, etc.) in the 
vicinity of the spray blocks. Under the sponsorship of the USDA Expanded 
Douglas-fir Tussock Moth Research and Development Program, work began in 
early 1977 on the development of a technical data base for use in applying 


mathematical dispersion models to USFS spray projects and the refinement 


and adaptation of existing models to predict spray behavior above and 
within forest canopies. This work (Dumbauld, Rafferty and Bjorklund, 
1977) resulted in the development of the CBG computerized spray dispersion 
model and the comparison of model predictions with measurements made 
during selected USFS spray programs. The CBG computer program has recently 
been applied in the development of optimum swath widths, application rates, 
and aircraft altitudes for the conduct of a pilot project in the Withla- 
coochee State Seed Orchard Project (Rafferty and Dumbauld, 1980) and is 
currently being used in evaluating the results of the project. The USFS 
(Dumbauld, Bowman and Rafferty, 1980) has also sponsored the use of the 
CBG model to determine optimum swath widths and spray altitudes for 
specific spray aircraft planned for use in the Maine 1980 Spray Program. 
A general discussion of spray behavior in and above canopies can be found 
in the report describing research funded by the USDA Expanded Douglas Fir 
Tussock Moth Research and Development Program (Brookes, Stark and Campbell, 
OTS) re 

The CBG computer based model is comprised of two major parts. 
The first part simulates the deposition of spray drops at the top of a 
forest canopy. The second part, which is based on a model developed by 
Grim and Barry (1975), simulates the penetration of spray drops into the 
forest canopy. This report describes a new computer program developed 
under the sponsorship of Methods Application Group, USFS (MAG) to extend 
the usefulness of the modeling concepts and applicability of the CBG 


computer based model. 


a2 PURPOSE 


The primary purpose for the development of the FSCBG computer 
program described in this report was to extend the CBG computer based model 
to include the effects of spray drop evaporation on the spray deposition 
patterns and drift concentration/dosage levels. Because the program 


changes for including evaporation effects were extensive, the decision was 


made to construct the new program in a modular fashion to ensure that 
further updates in the various calculation procedures could be more 
easily accomplished. Also, many calculations for deriving input para- 
meters, such as the gravitational settling velocity of drops, are per- 


formed automatically, on user option, by the new program. 


eS GENERAL DESCRIPTION OF THE FSCBG COMPUTER PROGRAM 


The FSCBG computer program implements an aircraft wake-settling 
velocity model, line-source dispersion models, a drop evaporation model 
and a canopy penetration model to calculate the concentration, dosage 
and deposition above and below a forest canopy downwind from multiple 
aircraft spray lines or swaths. The program is written in FORTRAN IV 
and is designed for implementation on a UNIVAC 1108 computer, although 
the program is easily adapted to other computer systems. The computer 
program requires 32K words of core storage. Figure 1-1 is a schematic 
diagram showing the major components of the FSCBG program. As indicated 
in the figure, program control parameters, as well as source and meteoro- 
logical data are input to the program via card image. The input require- 
ments of the program are discussed in detail in Section 3. Definitions 
of the meteorological and source inputs and their relationships to the 
various program models are given in Section 2. Figure 1-1 shows that 
the user may select a number of model options. The Above Canopy and 
Below Canopy Evaporation Model option allow the user to specify whether 
or not drop evaporation is to be included in a calculation of concentration, 
dosage and deposition. If the user selects a drop evaporation option, 
the program automatically calculates the change in drop diameter with 
time, using a polynomial expression fitted either to empirical data or 
theoretical calculations. The Wake-Settling Velocity Model option 
permits the user either to input the aircraft wake-settling velocity 
directly or to input the aircraft weight, wing span and ground speed 


which the program uses to calculate the aircraft wake-settling velocity. 
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Schematic diagram of the FSCBG computer program. 


The Canopy Penetration Model option permits the user to calculate the 
fraction of material presented at the canopy top, penetrating the canopy 
and reaching up to 10 levels within and below the canopy, including the 
ground. Finally the Dosage and Concentration Dispersion models permit 
the user to calculate dosage, concentration and area coverage of user 
specified dosage and concentration levels above the canopy and below the 
height of elevated inversions. The Deposition Dispersion Model allows 
the user to calculate deposition and deposition area-coverage at any 
height below the height of elevated inversions, including heights below 
the canopy top when the model is used in conjunction with the Canopy 


Penetration Model. In its present configuration, the FSCBG computer 


program can calculate the dosage, concentration and deposition at a 
maximum of 737 receptor coordinates downwind from a maximum of 100 line 


sources. 
Le& REPORT CONTENTS 


Section 2 of this report describes the aircraft wake model, 
dispersion models, evaporation model, canopy penetration model, as well as 
the meteorological and source inputs required by the models. User's 
instructions for the FSCBG program are contained in Section 3. An 
example problem illustrating use of the FSCBG program is given in Section 
4, Finally, a FORTRAN listing of the FSCBG program is contained in 


SecEi1on 5. 


SECTION 2 


DESCRIPTION OF THE FSCBG MODELS 
AND MODEL INPUT REQUIREMENTS 


Jacl. SPRAY DISPERSION MODELS 


The spray dispersion models incorporated within the FSCBG com- 
puter program are similar to the aircraft spray models described by 
Dumbauld, Rafferty and Cramer (1976). These models were developed pre- 
viously for the U. S. Army Dugway Proving Ground and are based on the 
generalized diffusion model concepts described by Cramer, et al. (1972). 
The models were also used in the CBG computer program and the mathematical 
basis for the models is given in Appendix B of the report by Dumbauld, 
Rafferty and Bjorklund (1977). In this report we describe the models 
used in the FSCBG program to calculate peak concentration, dosage and 
deposition downwind from nearly-instantaneous elevated line sources 
oriented at an arbitrary angle with respect to the mean wind direction. 
In the models, the axis of the spray cloud is assumed to be inclined 
from the horizontal plane by an angle that is proportional to Oe /u), 
where V, is the gravitational settling velocity for the j ch Patbe size 
we and u is the mean cloud transport speed. When evaporation is 
negligible, this inclination angle is invariant with distance from the 
line source. For evaporating drops, the angle changes with distance 


from the source because MM depends on the drop size. It is assumed that 


drops dispersed upwards by turbulence are totally reflected downwards at 
the top of the surface mixing layer, but the fraction of drops reflected 
at the ground is a variable input parameter for each ae category. The 
models use the Cartesian coordinate system shown in Figure 2-1 for a 


line source of length L at a height H. and a calculation point at R(e, 4, Zz). 


Dosage Models. The dosage (mass x time/unit volume) above 


the canopy is the sum of contributions from drops in the dispersing cloud 


FIGURE 2-1. Schematic plan view showing the line source geometry with 
respect to a calculation point at R(e, 6, z) for a wind 
direction 6. 


and from vapor produced by drop evaporation. The above canopy dosage 


for the contribution from drops is given by the expression 


where 
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The following parameters used in the preceding equations are 


based on meteorological measurements or inferred from meteorological 


measurements: 


On = standard deviation of the wind azimuth angle in 
radians 


alah 


k = constant relating On and OF 


= o'/o! (2-20) 
On = standard deviation of the wind elevation angle in 
radians 
Hi = depth of the surface mixing layer below a capping 
inversion 
u = mean transport wind speed above the canopy 
Q = angle between a line perpendicular to the line 
source and the mean wind direction (see Figure 2-1) 
Au = vertical wind-speed shear (see Equation (2-30)) 
The following parameters are source inputs required for use in the 
model: 
Q = total source strength emitted along the length 
L of the line source 
H = aircraft flight altitude above the ground 
V, = gravitational settling velocity for the median 
drop by mass in the qth drop-size category 
a = fraction of the total source strength in the 


3th drop-size category 


KZ 


V5 = reflection coefficient for the median drop by 
mass in the jth drop-size category 


95 = standard deviation of the cloud distribution at 
the distance Xp 
L = length of the line source 


As noted above, the vapor emitted from evaporating drops also con- 
tributes to the total dosage above the canopy. In the program construct, 
the evaporation model determines the mass of vapor emitted in Ah height 
units along the cloud axis as it descends to the canopy. The distance down- 
wind from the line source where the vapor is emitted is given by ut 
where t is the time from the evaporation model when the cloud axis passes 
through the midpoint of the height interval Ah. The source dimension of 


the vapor cloud emitted over the Ah height interval is 


(2-21) 


under the assumption that the vapor is rectangularly distributed at the 
point of emission. The vapor dosage is then calculated from Equation (2-1) 
with V; set equal to zero, V5 set to equal unity, and Xp set equal to ut. 
Finally, the total dosage at the calculations point is determined in the 


program by summing the contributions from the drops and vapor clouds. 


As inspection of Equation (2-1) and the definitions of the param- 
eters appearing in Equation (2-1) show, the model cannot be used when the 
wind direction is exactly parallel to the line source (9§=90 degrees) or 
when the wind direction is exactly perpendicular to the line source 


(8=0 degrees). For the case when §=90 degrees, the program uses the 


1S} 


exact solution by making the following substitutions in Equation (2-1) 


through (2-19): 


cot8@ = Q 
(n - 6 ae (2-22) | 
cos@¢ | 
sind = 1 


and defining x' as the downwind distance from the upwind edge of the line 


source and y as the crosswind distance from the line source. For the case 


when §=0 degrees, the program also uses the exact solution 


= En! Pa be i = 
dD, on a, => {Vertical Term} {Lateral Term} (2-23) 


The Vertical Term is 


oo . , (2i8 -H'+z+V.x/u : ae i 2iH +H'+z-V.x/u 2 
£. ft exp | = ee ea ee 
j j 2 os j Z oO 
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(2-24) 
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i r Or ; . 
where, for convenience, the quantity 0 is defined to be equal to unity. 
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The Lateral Term is 


3 erf ee a + erf ae _ (2-25) 
v2 Oy V2 Oy 


In Equation (2-24), x is the downwind distance from the line source and in 
Equation (2-25), y is the crosswind distance from the center of the line 
source. In the program, the standard deviation of the crosswind distribution 


O and of the vertical distribution oF are defined by the expressions 


= J - 
o Oy (x) + or (2-26) 
i] 
Ciie= ‘A (x) + oO (2-27) 
Z k fo) 


Peak Concentration Models 


The peak concentration (mass per unit volume) which occurs at the 


calculation point is determined from the expression 


Dy u 
OL = Parc (2-28) 
D8 


where 0, is the standard deviation of the cloud in the alongwind direction 


defined by 


; 1/2 
ou (#2) + a (2-29) 


IES) 


and 


L{x} = alongwind cloud length 
_ 0.6 Au (x: an sind | (2-30) 
u 2 
Au = vertical wind-speed shear in the layer containing 


the cloud 


The use of the peak concentration models requires that the meteorological 


parameter Au be specified in the program input list. 


Deposition Models 


The deposition, expressed in units of mass per unit area, at 
the point (€, 6, 0) downwind from line sources at an angle § with the 


mean wind direction is given by the expression 


2 


J fe 
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FE 
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(Equation (2-31) is continued) 
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(2-31) 
(Cont. ) 


2 
Ze et 
3 KDr exp( P) el gt/2 (oe ‘) ert (zt/2 (4 - x) 
Y; 3/2 z bd 25 a0 28 
4E 
V. Xy 
ae + ae (2-32) 
u 
oll V.Bk 5 
coals +(n -—*-) cots (2-33) 
On A cosé 
2 
2a 6 
k B +(n - ar | (2-34) 
Le | ers ; 
=| A (n 2 =) cote (2-35) 
- cos 
A u 


7, 
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Oo = cos@ 

A u 

Ded 5 ye 
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and the remaining parameters, except a and b, are identical to those defined 
for the dosage model. The definitions of a and b for the dosage model given 
by Equations (2-14) and (2-15) are reversed in the case of the deposition 


model, i.e., 


ese ey | (2-39) 


a(x" + x, - &sind) (2-40) 


The deposition model that applies when the wind direction is exactly parallel 
to the line source can be obtained from Equations (2-31) through (2-38) 

using the substitutions defined by Equation (2-22). Deposition for the 

case when the wind is exactly perpendicular to the line source is given by 


the expression 


oD 
fe OC=V) 
Deposition = pe be eee {Vertical Term} 
Ds 2T Oy o, (xtxy) 
j= (2-41) 


{Lateral Term} 
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where the Vertical Term is 


Cre 


en Wes ee HE aN, ox/e & 
+ > a 2 Hoe B= exp |- = eee (2-42) 
j De a 2 or 
i=l 
V.xy ee 2 Bee ch Va eer : 
+ y, | 2in +H! +— Soe  ———EEE 
j m j = 2 \ of 


As before, 0° is defined to be equal to unity and x is the downwind distance 


| from the source. The Lateral Term in Equation (2-41) is defined by Equation 
(2-25) where y is the crosswind distance from the center of the line 

source and - and Cle are defined, respectively, by Equations (2-26) and (2-27). 
(2-27). 


Ziad EVAPORATION MODEL 


The FSCBG computer program has a number of options allowing the 
user to account for the effects of the evaporation of drops in the calcula- 
tion of dosage, concentration and deposition. The first step in the calcula- 
tion of evaporation effects is the specification of the time-rate change 
of the drop diameter for the initial j drop-size categories. The program 
user has the option of selecting an automated procedure for the theoretical 
calculation of the time-rate change in drop diameter, or of supplying 


values of the constants A, B and C in the quadratic equation 


2 
Deee= AB © aCe 2- 
F 5 j j (2-43) 


specifying the change of drop diameter with time t after release. 


IY) 


The theoretical calculation of the time-rate change of drop 


diameter in the FSCBG program is based on the expression 


dD 


Dae = ae + at At (2-44) 
where 
De . = final diameter of the median drop in the Pee drop-size 
>J category after the time increment At 
Dindg ett : ' ; .th 
D, j = initial diameter of the median drop in the j drop- 


size category 


Frossling (see Fuchs, 1959, p. 44) defines the change in diameter of 


water drops due to evaporation by 


8 
aieD: (4x10) M, D_ p, (e -e_) 
—i - Oa MEMES! IR Boel Ree ET (2-45) 
t M D, Py (ee) Vv 
where 
My = molecular weight of evaporating vapor from the drop (g mois) 
Mo = mean molecular weight of the resulting vapor-air mixture in 
the transfer path which is approximated by that of air (M,) 
in the FSCBG program (g moll) 
D = molecular apa ele of the evaporating vapor in air at the 
by drop temperature (cm s~1) 
D. = drop diameter (um) 
S ; =o 
Ps = air density (g cm) 
E ete =3 
Pp = density of the liquid drop (g cm ~) 
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e = partial pressure of the evaporating vapor at the drop sur- 
S face (mb) 


e = partial pressure of the evaporating vapor at an infinite 
distance from the drop (mb) 


P = air pressure (mb) 


fe = ventilation factor (dimensionless) 


| 
| The air density in the model is calculated from the relationship 


e, = 4 (10) (2-46) 


where R* is the universal gas constant (8.31432 j moder ora) and Ty is 
| given in terms of the mixing ratio r,. According to Beard and Pruppacher 
| CIOVAD): 
1g 
. oe) 
Ce! pee) (2-47) 
where 


Ty = air temperature (°K) 
m= WOnezioT Ace =e.) (2-48) 


The diffusivity De of the vapor into air at the temperature of the surface 


of the drop TL depends on the liquid being considered. For water drops, 


the FSCBG program uses the expression (Pruppacher and Rasmussen, 1979) 
A 1.94 P | 
Dd. =o O20 T pe (2-49) 
fe) A 


yall 


where Le 1S (27 Seto Pies ER is, 1013625: mb and dies is determined from the 


relation 
tp ae ot Cae. 
T. = Th — re Ee ae 
k R Ty 
where 
L = latent heat of vaporization (cal cae) 
Dy 2 
= 597.3 T. 
r 
-4 
a = 0.107 +:°3..67 =. 10 TL 
€ p's -l1 -1lo,-l 
k = thermal conductivity (cal cm s KK) 
= Ka [2 = Clk? = 02 Iq) (_/P,)] 
cae 107>|5 69 - 0.017 (T.-T | 
d : ; r fe) 
-5 
kK. = 10 [3.78 - 0.020 (T,. - 7) | 
Ty = (I, + T,)/2 


The partial pressure SF is defined by the expression 


575.0466 +: 312382291 T; + 1.296028 re 


Ss 93.51611 - qT, 


Z2 


(2-50) 


(2-51) 


(2-52) 


(2-53) 


(2-53) 


(2-54) 


(2-55) 


(2-56) 


(2-57) 


(2-58) 


and 


B = 1 (2-59) 


The vapor pressure at e, is also obtained from Equation (2-57) when 


T. is replaced by T, in Equation (2-58) and B is replaced by (RH/100), where 


r A 
RH is the relative humidity in percent. To find the drop temperature, 
vapor pressure at the drop surface and diffusivity, Equations (2-49) 


through (2-57) are solved by iteration in the FSCBG program. 


The ventilation factor E for water is given by (Pruppacher and 


Rasmussen, 1979) 


iS vee 2 1/3 


0.78 + 0.308 Sc’/> Re 51.4 
fF = (2-60) 
1. o;paties se Rar! 1.4 


Ww 
1) 
JA 


Re 


Ww 
i) 
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where 


Re = Reynolds number 


sig Ds 10h, (a0 (2-61) 


Uy 


Sc = Schmidt number 


= (2=62) 
Vv Pa 


H, = absolute viscosity of air (g ene Sty) 
3/2 


-2 ie 
563427 x 1105) ( A (2-63) 


(T,+296.16) 296.16 


The user can, on option, ask the FSCBG program to theoretically 


determine the time-rate change of drop diameter for non-water drops. 
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However, additional information is required by the program. When this 
option is selected, the drop temperature and vapor pressure at the drop 


are automatically calculated from the expression (Picot, 1979) 


A exp (3 - c/7,,) Ser a (ea) ey - 1.) |, - A exp (8 - c/n.) (2-64) 
Vv 


i 
where 

A = P _/760 
fe) 

Ce = molal concentration of the liquid (mol ene) 

k = thermal conductivity of the gas mixture surrounding 
the drops at the drop surface temperature 

B,C = constants obtained from tables expressing variations 


of vapor pressure with temperature (see, for example, page 
D-140, Handbook of Chemistry, 58th Edition, published by 
Chemical Rubber Co.) 


Values of k, Dd. and L for non-water drops must also be specified by the 
program user. The value of Qo for non-water drops can probably be set 
equal to zero. When the parameters required by Equation (2-64) have been 
specified, the program uses Newton's iteration procedures to determine 

the drop temperature and vapor pressure and the time-rate change of 
diameter, using the following expression (Fuchs, 1959) for the ventilation 


coefficient 
Sh (2-65) 


where 


Sh 


Sherwood number 


parse Re) (2-66) 


and @ is set equal to 0.3 in the FSCBG program as susgested by: Picot, C1979) 


24 


Finally, for both water and non-water drops, the results obtained 
using the evaporation model for above canopy calculations are fitted with 
Equation (2-43) by least-squares over the time period required for the drop 
to evaporate to a diameter of 5 micrometers or the settling velocity to 
reach 0.02 miSia aendche yer is greater. For calculations of evaporation 
effects below the canopy, Equation (2-43) is fitted to the evaporation 
model results over the time period required for the drop to travel from the 


top of the canopy to the ground. 


The following parameters required by the evaporation model are 


based on meteorological measurements: 


Ty = air temperature 
Pa = air pressure 
RH = relative humidity (for water drops) 


The following parameters are source inputs required by the evap- 


oration model: 


Mo = molecular weight of evaporating vapor 
Dy. = density of the liquid for non-water drops 
Dis) CANOPY PENETRATION MODEL 


Grim and Barry (1975) developed a mathematical model to calculate 
the percentage material of a given drop-size category j which, after 
entering the forest at the top of the canopy, is retained at various levels 
within the canopy. The model is based on a Monte Carlo technique where 


a large number of drops in each size category are passed along a trajectory 


2S) 


through a simulated forest with trees assigned to equal areas according 
to the density (stems per acre) in the forest being simulated. The drop 
trajectory is a function of the gravitational settling velocity ue 
which varies along the trajectory when the evaporation model is used, 
and the mean wind speed at various levels within the canopy. As a drop 
proceeds along the trajectory, each tree is randomly displaced within 
the assigned area in the plane of the horizon and calculations are made 
to determine if the drop intersects the tree envelope and, if an intersec- 
tion occurs, whether the drop strikes a tree element. When the drop 
strikes a tree element, a tally is recorded for the height interval 
within the canopy where the "hit" occurs and for all greater height 
intervals. Drops proceed along the trajectory until a hit occurs or 
until the trajectory intersects the ground. After the specified total 
number of drops in the size category have passed along the trajectory, 
the tally number within each height interval is divided by the total 
number of drops to obtain the percentage of drops reaching the given 
height interval. The total number of drops passed along the trajectory 
required to achieve a stabilized solution (percentage penetration) is a 
function of the steepness of the trajectory, with more drops being 
required for size categories with large settling velocities. In the 
FSCBG program, 500 drops in each size category are passed along the 


trajectory. 


Figure 2-2 is a schematic diagram showing an example drop tra- 
jectory and forest construct. The drop trajectory in each of k height 


intervals is defined by the following piecewise linear function: 


gia = (2-67) 
0.25 (k-1) H, 
tan >, ho. 


26 


*4yso10J 
peqetnuts pue Aioqo.0fe1q doip oTduexe ue SuTMOYS WeAseTp OTJeWsYdS 


“e-¢ Wands 


727 


Vie e 0 (2-68) 


fae 0. 25KH | a Somer arm 0 i aes Smee Jom (2-69) 
where 
. = angle, measured in radians, defining the declination of 
the trajectory with respect to the plane of the horizon 
“4 Ve 
tan = =, p < 1.4 radians 
eek 
= (2-70) 
1.4 radians ; > > 1.4 radians 
H. = height of the forest canopy 
aoe = mean wind speed in the aes height interval within the 


canopy 


In Figure 2-2 the trajectory is shown as a linear function in each height 
interval. When the evaporation model is used, the trajectory is curvilinear 
in each height interval. The number of trees placed along the trajectory 


path in the simulated forest is given by the expression 


= 2-71 
Nip x, /Ax, ( ) 
where 
x, = maximum horizontal travel distance in meters of the 
drop within the forest canopy 
Ax, = alongwind tree spacing within the simulated forest 


= 63.615/vD, (2-72) 
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dD. = tree density in units of stems per acre 


As shown in Figure 2-3, each tree stem is given the following 


location along the trajectory: 


Re (N, - n) Ax, + (R -0.5) Ax, ser Madya Zr vekeliore ore N, (2-73) 

Ver (R' - 0.5) Ax, (2-74) 
where 

n = tree number 


and R and R' are uniform random numbers between 0 and 1. The possibility 
of the drop trajectory intersecting the tree envelope is determined by 
comparing the distance ye with the radius of the tree envelope (W,/2) at 
the height 4 where the trajectory passes through the distance Xe The 
tree widths W.. are specified by the program user at one-meter intervals 
and the program calculates the radius (Ww, /2) at A by linear interpolation. 


If ve is greater than (W./2), no intersection occurs and the computer 


program proceeds to the next tree along the trajectory. 


If the drop trajectory intersects the tree envelope, the 
program calculates the probability that the drop impacts on a tree 


element from the expression 


Ba = Tae PRPEN®) (2-75) 
where 
: ; ae th 
E. = impaction efficiency of the tree element for the j 
J drop-size category 
PRPEN = probability of penetration for the population of drops 


and for a horizontal trajectory through the tree 
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FIGURE 2-3. Schematic diagram showing the coordinate location of the 
"random tree". 
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t = path length correction factor for a non-horizontal 
trajectory 
Ea regi 
> Sap 
cos >, Wa 
at (2-76) 
H. H. 
ES ° > -— 
W epee 
™m ™m 
Ne = maximum width of the tree envelope 


In the FSCBG computer program, the user has the option of entering 
the impaction efficiency Be for each drop-size category or permitting the 
program to calculate ae from the following empirical relationship, recom- 


mended by Grim and Barry (1975) and attributed to Sell: 


2384) x 10° p? Vv 
ee ee 
3 J 
Ek, = 2-77 
j ( ) 
1 - J Opg weg EAL 
J 
where 
v = impaction velocity (ms —) 
- 2 | 1/2 
co + Ue (2-78) 
s = diameter of the element on which the drop impacts 


(cm) 
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A particular drop of the population of drops is assumed to intersect the 
tree element when the value of P, from Equation (2-75) is greater than a 
uniform random number R between O and 1. Each tree is divided into ten 
height-class intervals, and an intersection with a tree element is 
recorded as a "hit" in the class interval in which Z occurs and in 


every higher class interval. 


The process described above is repeated for every drop passed 
along the trajectory and the final percentage of material penetrating to a 
given height interval determined by dividing the number of recorded hits 
in the height interval by the total number of drops in all Hes size 
categories passed along the trajectory. When the evaporation model is 
used, the mass or number of drops reaching a given level in the canopy is 


further adjusted to account for losses due to evaporation. 


The inputs required by the canopy penetration model are: 


PRPEN = probability of penetration 
= P : th ; : Re, 
UU... = mean wind speed in the k” height interval within the 
‘i forest canopy 
D. = tree density in stems per acre 
H. = tree height in meters 
W = tree width at one-meter height intervals 


The computer program also permits the user to simulate a multi-storied 
canopy of up to three tree heights, with a different value of PRPEN 
possible for each tree or story height. The authors vary the values of 


PRPEN as a function of the four foliage types illustrated in Figure 2-4. 


32 


Foliage Type II Foliage Type III Foliage Type IV 


Foliage Type I 


PRPEN=0.13 


PRPEN=0. 62 PRPEN=0. 38 


PRPEN=0. 87 


MS eeaD SSR. way 


pecs aE NS ; 
Sil TEC g Pearly 
eg EEE a ist Sas 


Gane 


i 5G als é 
CIS 


63 


Qa 4 


Diagram illustrating four foliage types suggested for use in classifying forests 
(after Grim and Barry, 1975). 


FIGURE 2-4. 


The output from the canopy penetration model includes, in addition 
to the percentage of material in the see drop size category penetrating to 
each of ten levels within the canopy, the maximum horizontal travel distance 
of the drop within the forest canopy (x,)- When the deposition models 
described in Section 2.1 are used in conjunction with the canopy penetration 
model to calculate deposition within the canopy, the distance from the 
line source to the target is adjusted so that the deposition model calculates 


the deposition at the top of the forest canopy at a distance x, upwind from 


4 
the target. The deposition at the target is then determined by multiplying 
the deposition calculated from the deposition models in Section 2.1 by the 
percentage of material reaching the calculation height from the canopy 


penetration model. 


2.4 AIRCRAFT WAKE MODEL (WAKE SETTLING VELOCITY MODEL) 


During a short time period after spray is released from an 
aircraft, the vortices formed by the wings and propellers of fixed-wing 
aircraft and by helicopter rotors principally control the growth of the 
spray cloud. Except for the lateral translation of the vortex system by 
a crosswind, the vortices also control the position in space of the 
spray cloud and the amount of material which may be deposited on the 
underlying surface directly below the flight path. There are a number 
of sophisticated models presented in the literature describing conventional 
aircraft wake action (C. duP. Donaldson and A. J. Bilanin, 1975) and 
associated computer programs (Hirsh, 1978). Wickens (1977) has used 
numerical simulation models which neglect atmospheric diffusion to 
describe the effects of aircraft wake vortices on aerial spray deposition. 
Currently, the USFS is cosponsoring an effort with the National Aeronautics 
and Space Administration (NASA) to further investigate wake effects on 


spray deposition (Eckblad, 1980). 


The relatively complex wake model programs currently available 


would, if incorporated with the dispersion evaporation and canopy penetration 
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codes in the FSCBG program, greatly increase the computation time required 
to estimate deposition patterns. From discussions with USFS personnel and 
others involved in the USFS-NASA program on wake effects, it appeared that 
simpler wake models and computer codes may soon be available. For these 
reasons, the FSCBG program has been designed so that wake computer codes 


can easily be modified and/or added to the program. Since the importance 


of wake effects was recognized, the version of the FSCBG code described 
i in this report contains an interim wake-effects model we have found 


useful in previous work (Dumbauld, Rafferty and Bjorklund, 1977). 


Prandtl and Tietjens (1934) express the sink rate W (m a) of 


a vortex system as 


| 8g W s 
| 0 (10 :) (2-79) 
| 3 2 
Teese Bein. 
| where 
| Ws = weight of the aircraft (kg) 
g = accelleration due to gravity (9.8 m eas 
Ons air density (g an 
b = aircraft wing span (m) 
Ve = aircraft speed (m Sis) 


Equation (2-79) is strictly applicable to fixed-wing aircraft. However, 
as pointed out by the Bell Helicopter Company (1966 and updates), helicopter- 


produced vortices resemble fixed-wing vortices at high forward speeds. 
In the spray dispersion models described in Section 2.1 above, 


the source parameters correspond to the spray cloud properties when the 


cloud has reached approximate atmospheric equilibrium. The vortex sink rate 
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calculated from Equation (2-79) and the observations that the vortex tube 
stops descending at a distance (b/2) above the surface (Jones, 1970) or 
above the forest canopy are used to calculate source parameters for the 
models. Figure 2-5 is a schematic diagram showing the geometrical 
considerations used in specifying the distance downwind Xp at which cloud 
stabilization occurs and the effective release height Hi - As shown 

in Figure 2-5, the cloud descent from the aircraft at height H forms the 
angle eanc- (w/u) with the horizontal when the gravitational settling 
velocity for the ge drop-size category at height Be is less than 

W. In general, the FSCBG program calculates the Hileeanee downwind from 
the flight path where the cloud reaches a height b/2 above the surface 


from the expression 


= = U _ 
Xp =u 2 (2-80) 
where 
t! = effective time for the cloud centroid to reach the 
J height b/2 
J = ° > 
(CH= (H/2)e= Ui ynGe Va se aah erin 
= (2-81) 
: < 
“45b/2 OSG SH 
: : , : - tn : 
Vig = gravitational settling velocity for the j drop-size 
Js category at the aircraft flight altitude 
: , .th 
CF .b/2 time for the drop in the j category to reach the 


height b/2 above the ground or canopy, calculated from 
the evaporation model 


As noted in Section 2.1, the dispersion models are formulated under the 


assumption that the cloud axis descends at an angle tan” * (V,/a) with the 
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horizontal. We have therefore assumed that the effective declination of 
the cloud axis downwind from an effective source height et is given by the 
is the gravitational settling velocity 


angle tant /u), where v, 


V5 3b/2 3b/2 
calculated by the program for the drop at the height b/2. Based on this 


assumption, the effective release height H! is 


Hed =) Hie 6b /2) e+ EV 2-82 
; if (b/2) ; ( ) 


tsb /2 


All gravitational settling velocities are calculated in the FSCBG program 
using a technique suggested by McDonald (1960). The graph given by McDonald 
expressing the relationship between the drag coefficient of spheres and 
their Reynolds number, has been fitted for Reynolds numbers less than 


2x10? and is used in the FSCBG program. 


PAS ADJUSTMENTS OF DEPOSITION CALCULATIONS FOR EVAPORATION 


When no evaporation occurs, the calculation of. deposition using 
Equations (2-31) through (2-32) is relatively straightforward because the 
gravitational settling velocity is invariant with time. However, when 
evaporation occurs both drop size and gravitational velocity vary with 
travel distance from the source and the deposition calculations become 
more complicated. To illustrate the procedures used by the computer program 
tc calculate deposition with evaporation occurring, assume that we wish 
to calculate prennd-teyel. deposition for the a drop-size category at a 
distance of 75 meters downwind from an aircraft release at a height of 
18.1 meters above the ground. The top of the canopy is at a height of 
12 meters above the ground. The program follows the stepwise procedure 


outlined below. 


(1) The program uses the above-canopy evaporation model de- 


scribed in Section 2.2 to establish a curve relating the 
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(2) 


(3) 


(4) 


gravitational settling velocity (V.) of the median drop 
in the gel size category as a ee of travel distance 
(x) from the source for above-canopy meteorology. A 
typical curve (for drop-size category 9 of the example 
problem explained in Section 4 and a relative humidity of 


50 percent) is shown in Figure 2-6. 


The above-canopy evaporation model is also used to cal- 
culate the height trajectory. of the drop as a function of 
travel distance. The curve for drop-size category 9 in 
the example problem for 50-percent humidity is shown in 
Figure 2-7. Note that the curve extends to negative 
heights because a drop in this category would have fallen 
this distance except for the effects of vertical disper- 


sion above the canopy. 


The below-canopy evaporation model is used to calculate 
the trajectory of the drop in the ce size category, 
using below-canopy meteorology and assuming that the drop 
entered the canopy with the settling velocity ie given 
by the above-canopy evaporation model at the canopy 
height. In addition to calculating the mass loss within 
the canopy, the below-canopy evaporation model calculates 
the travel distance (DISTM) of the drop from the time it 


enters the canopy until it strikes the ground. 


The program internally develops a curve relating the dif- 
ference between the receptor distance from the source and 
the travel distance in the canopy Corus for all the 
drop-size categories versus the settling velocity at the 
canopy top vie for each flan category. The curve for our 


example is shown by the nearly vertical line in Figure 
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FIGURE 2-6. Gravitational settling velocity V. versus distance from 


the source (see Section Zee 
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2-6, where the quantity CEST ED) - is plotted as X and Mie 


is plotted as rhe 


(5) The intersection of the curve developed in step (4) above 
and the above-canopy evaporation curve for drop-size cate- 
gory 9 shown in Figure 2-6 defines the distance from the 
source to the point where a drop in size category 9 would 
have to enter the canopy in order to reach the vicinity 
of our receptor at a distance of 100 meters from the source. 
Thus, for our example, the drop in size category 9 would 
enter the canopy at a distance of 87 meters downwind 


from the source. 


(6) The effective height H" of the cloud centroid for size 
category 9 at a distance of 87 meters from the source is 
obtained from the height versus distance curve shown in 
Figure 2-7. For our example, Figure 2-7 shows that H" 


is about -/4 meters. 


(7) The deposition model equations yield the corrected depo- 
sition values when an effective value of - is used in 
solving Equation (2-31) through Equation (2-42). Thus, 
the effective value of the gravitational settling velocity 
is 

u(—H! - H") 


x 


where H, is given by Equation (2-82) and, for our example, 


x is 87 meters and H" is -74 meters. 


The calculations then proceed in the usual fashion except for appropriate 


adjustments for mass loss due to evaporation. 
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Height of drop in size category 9 versus distance from the 


Source (see Section 2.5). 
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SECTION 3 


USERS INSTRUCTIONS FOR THE FSCBG 
COMPUTER PROGRAM 


Bleyll SUMMARY OF PROGRAM MODEL OPTIONS, DATA REQUIREMENTS AND OUTPUT 


Bilal Summary of FSCBG Program Model Options 


The program model options of the FSCBG computer program consists 


of four general categories: 


© Wake-Settling Velocity Model 


e Drop Evaporation Model 
e Canopy Penetration Model 
e Dispersion Models 


Each category is discussed separately below. 


Wake-Settling Velocity Model. The FSCBG computer program provides 


the option of allowing the user to input the wake-settling velocity or 
allowing the FSCGB program to calculate the wake-settling velocity, pro- 
viding the user inputs the weight of the aircraft, the ground speed of 


the aircraft and the wing span of the aircraft or helicopter rotor diameter. 


Drop Evaporation Model. The FSCBG program requires the user to 
input the spray drop-size distribution. The other program models use these 
data with or without drop-size changes due to evaporation. If the user 
chooses to use the drop-size distribution without changes due to evaporation, 
the drop size as well as all related drop information (settling velocity, 
etc.) is assumed to be constant with time. However, if the user chooses 
to apply evaporation to the drop distribution, the program will determine 
the rate of change of the drop size as well as all related drop information 


(settling velocity, etc.) as a function of time. If desired, the FSCBG 
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program permits the user to specify the rate of change of drop size for 
each drop-size category as a function of time. If this option is chosen, 
the program requires the user to specify, for each drop-size category 
above the canopy, the coefficients a, b, and c of the equation [ drop = 
atop. Er "at a » where drop is the drop diameter in micrometers 

and t is the time in seconds from release. Also, the program permits 

the user to specify a different set of coefficients (a, b and c) for each 
drop-size category giving the time rate of change of drop size within 


the forest canopy. 


Canopy Penetration Model. The user may run the FSCBG program 
with or without a forest canopy. A forest canopy may be any vegetation 
that extends above the ground, where the canopy can be described using up 
to 3 vegetative height classes. The canopy penetration model calculates 
the percentage of material in each drop-size category that reaches each 
tenth of the canopy height. Also, the canopy penetration model determines 
the horizontal path length of the trajectory of each drop-size category. 
This information can be optionally printed and is supplied to the dispersion 
model (if used) for application in the calculation of dosage, concentration 


and/or deposition. 


Dispersion Models. The user has the option of using the dosage, 
concentration or deposition models or any combination of these models. The 
program on option will calculate the dosage, concentration and/or deposition 
at the specified receptor points and at the specified height. At present, 
dosage and concentration calculations are restricted to heights at or above 
the canopy height. There are no similar height restrictions on deposition 


calculations. 
oude2 Data Input Requirements 


This section provides a description of all input data parameters 
required by the FSCBG program. The user should note that some input param- 
eters are not read or are ignored by the program depending on the values 


control parameters have been assigned by the user. 
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Program Control Parameter Input Data. These parameters provide 


user-control of all program options. 


Parameter 
Name 


TITLE Titling Information -- This parameter is an array of a maxi- 
mum of 80 characters specifying the case title. If not 


input, the program defaults to a blank title. 


ISW(1) Wake-Settling Velocity Option -- This parameter specifies 
whether the wake-settling velocity is to be input by the 
user or is to be calculated by the program. If this param- 
eter is set equal to "0" or omitted, the FSCGB program 
assumes the user desires to input the wake-settling veloc- 
ity via the parameter WAKVEL below. If this parameter is 
set equal to "1", the program calculates the wake-settling 
velocity from the input parameters ARCRWT, WNGSPN, ARCRSP 
and AIRDEN, defined below. 


ISW(2) Evaporation Model Option -- This parameter specifies whether 
or not evaporation is to be applied to the drop distribution 
or whether the user is to input the rate of evaporation of 
the drops. If this parameter is set equal to "0" or omitted, 
the program assumes the drops do not evaporate, and passes 
the drop-size distribution to subsequent program phases 
unchanged. Also, all equations related to drop size and time 
assume the drop size is constant with time. If this param- 
eter is set equal to a value of "1", the program calculates 
the evaporation rate for each drop-size category. The pro- 
gram automatically calculates the drop size as a function of 
time and the amount of material converted to vapor. The 
evaporation model is applied to the drop distribution 


both above and below a forest canopy and the evaporated 
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Parameter 
Name 


TSW(2) 
(Cont. ) 


ISW(3) 


ISW(4) 


drop distributions are passed to subsequent program phases. 
If this parameter is set equal to "2", the FSCBG program 
assumes the user desires to input the coefficients a, b 

and c of the equation giving the drop diameter as a 
£unction-of time, e.g. 5 D-= a+b ¢ tit.c Ee. One set 

of the coefficients a, b and c is required for each drop- 
size category for both above and below a forest canopy. The 
variable D in the above equation is the drop diameter in 


micrometers and t is the time after release. 


Canopy Penetration Model Option -- This parameter controls 
execution of the canopy penetration model. If this parameter 
is set equal to "0" or omitted, a forest canopy is not con- 
sidered in the FSCBG calculations. All input parameters re- 
lating to a forest canopy are ignored by the program. If 
this parameter is set equal to "1", the program assumes there 
is a forest canopy. The forest canopy may be an actual forest 
or may be any above-ground vegetation. The canopy penetra- 
tion model determines the percentage of the drop material 
from each size category that reaches 10 height levels 

within the forest canopy, including the ground. If this 
parameter is set equal to "2", the canopy penetration model 


calculations are printed as well as calculated. 


Dispersion Models Option -- This parameter controls execu- 
tion of the dispersion models. If this parameter is set 
equal to "0" or omitted, the dispersion models (dosage, 
concentration, deposition) are not executed. All dispersion 
model related parameters are ignored by the program. If 


this parameter is set equal to "1", the dispersion models 
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Parameter 
Name 


ISW(4) 
(Cont. ) 


ISW(5) 


ISwW(6) 


ISW(7) 


ISW(8) 


given by ISW(6), ISW(7) and ISW(8) are executed and the 


calculations printed. 


Print Evaporation Model (ISW(2)) Calculations Option -- This 
parameter controls printing of calculations selected under 
Isw(2). If this parameter is set equal to "0" or omitted, 
calculations under ISW(2) are not printed. If this param- 


eter is set equal to "1", calculations under ISW(2) are 


_ printed. 


Dosage Model Option -- This parameter controls dosage calcu- 
lations. If this parameter is set equal to "0" or omitted, 
dosage is not calculated. If this parameter is set equal to 
"1", dosage at the height specified by Zz is calculated and 
printed. Note that dosage calculations can only be made at 
heights Z from the top of the forest canopy to the mixing 
layer height HM. 


Concentration Model Option -- This parameter controls concen- 
tration calculations. If this parameter is set equal to "0" 
or omitted, concentration is not calculated. If this param- 
eter is set equal to "1", concentrations at the height speci- 
fied by Z are calculated and printed. Note that concentra- 
tion calculations can only be made at heights Z from the 


top of the forest canopy to the mixing layer height HM. 


Deposition Model Option -- This parameter controls deposi- 
tion calculations. If this parameter is set equal to "0" 

or omitted, deposition is not calculated. If this parameter 
is set equal to "1", deposition at the height specified by 


Z is calculated and printed. There are no restrictions 
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Parameter 
Name 


ISW(8) 
(Cont. 


ISW(9) 


ISW(10) 


ISW(11) 


on the calculation height Z for deposition calculations, 


e.g., deposition can be calculated within the forest canopy. 


Dosage, Concentration and Deposition Output Mass Units -- 
This option assumes the user has input the release rate Q 
in units consistent with the definitions of Q and with the 
parameter SWATH below. All units conversions are handled 


automatically by the program. The possible output units 


ares 
TSW(9) Output Units 
O or omitted drops 

1 micrograms 

7 milligrams 

S) grams 

4 ounces 

5 pounds 
Dosage Output Time Units Option -- This option assumes the 


user has input all parameters dependent on time in units of 


seconds. The possible output units are: 


ISW(10) Output Units 
O or omitted seconds 
1 minutes 


Dosage, Concentration and Deposition Output Volume or Area 
Length Units Option -- This option assumes the user has in- 
put all parameters dependent on length in units of meters. 


The possible output units are: 
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Parameter 


Name 
ISW(11) ISw(11) Dosage Concentration Deposition Area-Coverage 
(Cont. ) 
0 or omitted ae a ma ae 
1 ft? ft? et” fe 
3 3 
2 T™. ™ acre acre 
| 3 ft? ge acre acre 
| 4 a te hectare hectare 
| Isw(12) Area-Coverage Option -- This option enables the user to cal- 
culate the total horizontal area coverage for up to 10 levels 
of dosage, concentration and/or deposition. If this param- 
eter is set equal to "0" or omitted, area-coverage is calcu- 
lated for dosage, concentration and/or deposition, depending 
js on ISW(6), ISW(7) and ISW(8), respectively. 
ISW(13) - Reserved for Future Options. 
ISW(20 


Necessary FSCBG Program Input Parameters. These parameters are 


used by more than one of the FSCBG program models. 


Parameter 
Name 


IFWATR Spray Liquid Base -- This parameter is a switch used to 
inform the FSCBG program when the theoretical calculation 
for water drop evaporation is to be used. If this parameter 
is set equal to "1" or omitted, the program assumes the theo- 
retical drop evaporation equations for water are to be used. 
If this parameter is set equal to "2", the program assumes 
the spray liquid is not water, but calculates a theoretical 


evaporation rate requiring these additional input parameters - 
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Parameter 
Name 


TFWATR 
(Cont. ) 


WNGSPN 


WSOCAN 


HGTCFT 


DENLIQ 


AIRTPO 


AIRTPU 


ATRDEN 


DFUSIV, HETLAT, CONMOL, THERMC, VAPINF, BCONST and CONST 


defined below. 


Aircraft Wing Span -- This parameter specifies the aircraft 


wing span or helicopter rotor diameter in meters (m). 


Wind Speed -- This parameter specifies the wind speed in 


meters per second (m ona) above the canopy. 


Height of Aircraft -- This parameter specifies the height 


above ground in meters (m) of the spray aircraft. 


Density of the Drop Liquid -- This parameter specifies the 
density of the spray drop liquid in grams per cubic centi- 
meter (g cm >), If this parameter is omitted from the input 
data, the program defaults to a density of 1.0 gram per cubic 


centimeter. 


Air Temperature Above the Canopy -- This parameter specifies 


the air temperature in degrees Celsius Cc) above the canopy. 


Air Temperature Below the Canopy -- This parameter specifies 
the air temperature in degrees Celsius (°C) below the canopy. 
If this parameter is omitted from the input data, the program 


defaults to the value input to AIRTPO. 


Air Density -- This parameter specifies the air density in 
grams per cubic centimeter (g cm). If this parameter is 
omitted from the input data, the parameters AIRMOL, AIRPRS 
and VAPINF defined below, must be input to the FSCBG program. 


50 


Wake-Settling Velocity Model Input Data. These parameters are 


used in calculating the aircraft wake-settling velocity. The parameter 
ISW(1) specifies which of the parameters described in this section are 


required by the program. 


Parameter 
Name 
Parameters Required When (ISW(1) Equals "0" or is Omitted 

WAKVEL Wake Settling Velocity -- This parameter specifies the wake- 
settling velocity in meters per second (m aie )a 
Parameters Required When ISW(1) Equals "1" 

ARCRWT Aircraft Weight -- This parameter specifies the weight of 
the aircraft in kilograms (kg) for the calculation of the 
wake-settling velocity. 

ARCRSP Aircraft Ground Speed -- This parameter specifies the ground 


speed in meters per second (m ae) of the spray aircraft for 


the calculation of the wake-settling velocity. 


Evaporation Model Input Data. These data are parameters used in 


calculations within the Evaporation Model. The parameter ISW(2) specifies 


which of the parameters in this section are required by the program. 


Parameter 
Name 
Parameters Required When ISW(2) Equals "0", "1" or is Omitted 
DRPUPR Upper Limits for the Drop-Size Categories -- This parameter 


is an array specifying the upper limit of each drop-size 


ill 


Parameter 
Name 


DRPUPR 
(Gont=) 


DRPLWR 


ATRPRS 


AITRMOL 


VAPMOL 


category in micrometers (um) for up to 20 drop-size cate- 
gories. These values must be input in descending order 

of diameter and the number of drop-size categories is deter- 
mined by the program from the number of non-zero input 
values. If ISW(2) is equal to "0" or omitted, this array 


may be the mean diameter. 


Lower Limits to the Drop-Size Categories -- This parameter 
is an array specifying the lower limit of each drop-size 
category in micrometers (um) for the same number of drop- 
size categories input to the array DRPUPR. The lower limit 
of any drop size category cannot be zero ("0"). If the 
mean diameter is input to the array DRPUPR, this array is 


omitted from the input data. 


Parameters Required When ISW(2) Equals "1" 


Air Pressure -- This parameter specifies the barometric 
pressure in millibars (mb) at the site altitude. If 

this parameter is omitted from the input data, the program 
uses 1013.25 millibars as a default value. If the param-— 


eter AIRDEN has been omitted, this parameter is required. 


Molecular Weight of Air -- This parameter specifies the 

molecular weight of air. If this parameter is omitted from 
the input data, the program uses 28.9644 as a default value. 
If the parameter AIRDEN has been omitted, this parameter is 


required. 


Molecular Weight of the Vapor -- This parameter specifies 


the molecular weight of the vapor from the evaporating 


az 


Parameter 
Name 


VAPMOL 
«Cont*)) 


RELHMO 


RELHMU 


DFUSIV 


HETLAT 


drops. If this parameter is omitted from the input data, the 


program uses 18.015 as a default value. 


Relative Humidity Above the Canopy -- This parameter speci- 
fies the relative humidity above the canopy in percent (%). 
This parameter is required if AIRDEN equals zero or is omit- 


ted and VAPINF is omitted from the input data. 


Relative Humidity Below the Canopy -- This parameter speci- 
fies the relative humidity below the canopy in percent (%). 
If this parameter is omitted from the input data, the pro- 
gram uses the humidity above the canopy (RELHMO) as a default 
value. This parameter is required if AIRDEN equals zero or 


is omitted and VAPINF is omitted and a canopy is being used. 


Parameters Required When ISW(2) Equals "1" and IFWATR Equals "2" 
Parameters nKequired wnen tow (4) Equals 41) and tfWwAtn EqualS 4 


Diffusivity of Evaporating Vapor -- This parameter is required 
for liquids other than water (IFWATR="2"') and specifies the 
diffusivity of the evaporating vapor into the air at the drop 
temperature in square centimeters per second Gone ie liek 
this parameter is omitted from the input data and IFWATR 
equals '2", the program calculates DFUSIV, assuming the 


liquid is similar to water via the equation 


DFUSIV = .211- ((1,+273.16)/273.16)"*7*+(1013.25/ATRPRS) (3-1) 


where Ty is the drop temperature approximated by the program. 


Latent Heat of Vaporization -- This parameter is required 


only for liquids other than water (IFWATR="2") and specifies 
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Parameter 
Name 


HETLAT the latent heat of vaporization at the drop temperature in 
icone.) calories per mole (cal mole‘), If this parameter is omitted 
from the input data and IFWATR equals "2", the program cal- 
culates HETLAT, assuming the liquid is similar to water, 


via the equation 
HETLAT = 597.3+ (273.16/ (1,273.16) )*+VAPMOL (352) 


where: 


A 0.107+3.67x10" "+ (T,+273.16) 


TD 


drop temperature (°C) approximated by the program 


CONMOL Molal Concentration -- This parameter is required only for 
liquids other than water (IFWATR="2") and specifies the 
molal concentration of the air-liquid mixture in moles per 
cubic centimeter (moles cm). If this parameter is omitted 
from the input data and IFWATR equals "2", the program calcu- 


lates CONMOL, assuming the liquid is similar to water, via 


the equation 


CONMOL = AIRDEN/AIRMOL (3-3) 


THERMC Thermal Conductivity -- This parameter is required only for 
liquids other than water (IFWATR="2") and specifies the 
thermal conductivity of the vapor into air at the drop tenmn- 
perature in calories per second centimeter degree Kelvin 


si 


= -l1 o,-L . Sea ET 
(cal s cm K “). If this parameter is omitted from 


the input data and IFWATR equals "2", the program calculates 


THERMC, assuming the liquid is similar to water, via the 


equation 
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Parameter 
Name 


THERMC 
(Cont. ) 


VAPINF 


BCONST 
& 
CCONST 


THERMC = a(1- (1.17-1.02 (5/a))) VAPINF/ALRPRS (3-4) 
where 
-5 -7 
AS= 25 69x Oee 17x10) DRE IMP 
B = 3.78x10 7-+2.0x10~’+DRPTMP 
DRPTMP = calculated drop temperature (°K) 
VAPINF = vapor pressure of vapor at infinity (mb) 
AIRPRS = barometric pressure (mb) 
Vapor Pressure at Infinity -- This parameter is required for 


liquids other than water (IFWATR="2"') and is also required if 
AIRDEN equals zero or has been omitted from the data deck. 

This parameter specifies the vapor pressure of the evaporating. 
vapor at infinity in millibars (mb). If VAPINF is input 
greater than or equal to zero, this value is used. However, 
if this parameter is omitted from the input data, the pro- 
gram calculates the value of VAPINF from the air temperature 


(AIRTPO,AIRTPU) and relative humidity (RELHMO,RELHMU) . 


Vapor Pressure Equation Constants -- These two parameters are 
constants used in the equation that describes the vapor pres- 
sure as a function of temperature of the non-water (IFWATR="2") 


liquid via the expression 
vapor pressure (in Hg) = exp (BCONST-CCONST/T,,) (3-5) 
where Ty is the drop temperature (Ge) approximated by the 


program. The default values for these parameters are 21.07 


and 5249.9 for BCONST and CCONST, respectively. 
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Parameter 
Name 


Parameters Required When ISW(2) Equals "2" (User Specifies 
Time Rate of Evaporation). 


DAU, DBU Time Rate of Change of Drop Size Above Canopy -- These three 
arrays specify the coefficients of the equation describing 


drop size as a function of time above the canopy 
D(micrometers) = DAU(J)+DBU(J) *T+DCU(J)°T? (3-6) 


where D is the diameter of the drop in micrometers (um), T 

is the time in seconds after the spray is released and J is 
the index over the drop-size categories. There are a maximum 
of 20 values possible for each array. The number of drop- 
size categories is determined by the program by the number 

of input values of DAU. The values in each array must be 


ordered in descending order of drop-size category. 


DAL, DBL Time Rate of Change of Drop Size Below Canopy -- These three 
arrays specify the coefficients of the equation describing 


drop-size as a function of time below the canopy 
D(micrometers) = DAL(J)+DBL(J) *T+DCL(J) °T? (3-7) 


where D is the diameter of the drop in micrometers (um), 

T is the time in seconds after the spray is released and J 
is the index over the drop-size categories. There are a 
maximum of 20 values possible for each array. The values 
in each array must be ordered in descending order of drop- 
size category. If any values for these three arrays are 
omitted from the input data, the program defaults the re- 


pective value from the arrays DAU, DBU and DCU. 
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Canopy Penetration Model Input Data. These parameters are used 


in calculations performed by the Canopy Penetration Model and are required 


by the program only if the parameter ISW(3) is set equal to "1". 


Parameter 
Name 


WSIN44 Wind Speed in Top Quarter of Canopy -- This parameter speci- 
fies the wind speed in meters per second (m oa) for the top 
quarter of the canopy. If this parameter is omitted from the 
input data, the program uses the wind speed above the canopy 


(WSOCAN) for WSIN44. 


WSIN34 Wind speed in the Third Quarter of Canopy -- This parameter 
specifies the wind speed in meters per second (m en) for the 
third quarter of the canopy. If this parameter is omitted 


from the input data, the program uses the value of WSIN44. 


WSIN24 Wind Speed in the Second Quarter of Canopy -- This parameter 
specifies the wind speed in meters per second (m a) for the 
second quarter of the canopy. If this parameter is omitted 


from the input data, the program uses the value of WSIN34. 


WSIN14 Wind Speed in the Bottom Quarter of Canopy -- This parameter 
specifies the wind speed in meters per second (m BU) for 
the bottom quarter of the canopy. If this parameter is omit- 
ted from the input data, the program uses the value of 


WSIN24. 


HGTCAN Height of Canopy -- This parameter is an array specifying 
the height in whole meters (m) of up to three canopy (vegeta- 
tion) height categories (multi-storied canopy). These 


values must be input in descending order of height. The 


Sy} 


Parameter 
Name 


HGTCAN number of canopy height categories. is determined by the 
one program from the number of non-zero input values of HGTCAN. 


There are a maximum of three height categories. 


PRBPEN Probability of Drop Penetration -- This parameter defines 
an array specifying the probability of drop penetration for 
a horizontal trajectory through the tree (stem) for up to 
three canopy height categories. These values are input as 
fractions in the range of "0" to "1", depending on the den- 
sity of the vegetation and in the same order as values input 


to the array HGTCAN. 


TREDEN Tree Density -- This parameter is an array specifying the 
density of trees per acre (stems acre +) for each canopy height 


category defined by HGTCAN. 


COLEFF Collection Efficiencies -- This parameter is an array that 
either specifies the collection efficiencies for each drop- 
size category or specifies the size of the vegetative ele- 
ments for each canopy height category, depending on whether 
the values are negative or positive. If the values input 
to this array are positive, the program assumes they are the 
collection efficiencies for each drop-size category. If the 
values input to the array COLEFF are negative, the program 
assumes the absolute values are the diameters in centimeters 
of the vegetative element on which the drop impacts for each 
canopy-height category. The collection efficiencies are then 
internally calculated using the following empirical relation- 
ship, recommended by Grim and Barry (1975) and attributed to 
Sell 
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Parameter 


Name 
COLEFF 
Keone ey 2.8x10-*+D,+u/s ; EFR SS <1 
EFF = (3-8) 
it : BER 
where 
u = impaction velocity in meters per second (m sly 
D, = drop diameter for the jth drop-size category 
J in micrometers (Um) 
s = diameter in centimeters of the element on which the 
drop impacts (input as a negative value in COLEFF) 
If this array is omitted from the input data, the program 
will default to a vegetative element size of 13 centimeters 
by setting the array equal to all values of -13. 
TREENV Tree Envelope -- This parameter is an array specifying the 


width of the tree height class in meters (m) at each meter 
of tree height. Values are input in ascending order of the 
height of the tree class, beginning at one meter. There are 
exactly HGTCAN(J) values input for each canopy class J. The 
first HGTCAN(1) values are for height category 1, the second 
HGTCAN(2) values are for height category 2 and the last 
HGTCAN(3) values are for height category 3. 


Dispersion Model Input Data. These parameters are used in the 


calculation of dosage, concentration and/or deposition and are required by 


the program only if ISW(4) is equal to "1". 
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Parameter 
Name 


NSOURC 


NXPNTS 


NYPNTS 


NXYPNT 


Number of Line Sources -- This parameter specifies the 
number of line sources or aircraft spray lines to be 
used. The program is capable of processing a maximum of 
100 line sources. If this value is input as "0" or 
omitted from the input data, the program defaults to a 


value of "1". 


Number of X (East-West) Receptor Coordinates -- This 
parameter specifies the number of X (east-west) receptor 
coordinates in the rectangular receptor grid system. The 
maximum number of X receptors depends on the values of 


NYPNTS and NXYPNT and is limited by 


O < NXPNTS+NYPNTS+2*NXYPNT < 100 and 
and (3-9) 
O < NXPNTS*4NYPNTS+NXYPNT S81) 


Number of Y (North-South) Receptor Coordinates -- This param- 
parameter specifies the number of Y, (north-south) receptor 
coordinates in the rectangular receptor grid system. The max- 


imum number of Y receptors is given by Equation (3-9). 


Number of X,Y Discrete (Arbitrarily Placed) Receptor 
Coordinates -- This parameter specifies the number of X,Y 
discrete receptor coordinates. The maximum number of 


discrete coordinates is given by Equation (3-9). 


Emission Rate -- This parameter specifies the spray 
emission rate for the line sources in units of grams per 
meter, or gallons per acre, depending on the input param- 


eter SWATH. If the parameter SWATH is greater than zero, 
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Parameter 
Name 


Q 
(Cont. ) 


SWATH 


the program assumes Q is in gallons per acre and the area 
sprayed is sprayed in lines of an equal distance (SWATH) 
apart. If the spray lines are not uniform or the area is 
not regular, input Q in units of grams per meter and set 
the parameter SWATH equal to "0" or omit SWATH from the 
input data. If the parameter Q is input as "0" or omitted 
from the input data, the program defaults to a value of 
"1", The equation used to convert gallons per acre to 


grams per meter is given by 


Q(g/m) = £9352938(acre+cn /gaillem-) -DENLTO(s/em >): 
SWATH (m) °Q(gal/acre) 


where 
DENLIQ = input liquid density 
SWATH = input swath width 
Q = input emission rate (gal/acre) 
Distance Between Spray Lines -- This parameter specifies 


the distance between uniform spray lines in meters. If 
the input of units of Q are in gallons per acre, then 
SWATH must be input as a value greater than zero. However, 
if the input units of Q are grams per meter then SWATH 


must be input as "0" or omitted from the input data. 


Height of Calculations -- This parameter specifies the 


height above ground in meters at which the calculations 


of dosage, concentration and deposition are to be performed. 


The FSCBG program at present does not have the capability 


to calculate dosage or concentration below the canopy 
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Parameter 
Name 


(Cont. ) 


SIGAP 


SIGEP 


TAU 


TAUO 


height. However, deposition can be calculated at any 
height below the mixing depth HM. If this parameter is 
set to a value of "0" or omitted from the input data, the 


program assumes ground-level calculations are to be made. 


Standard Deviation of the Wind Azimuth Angle -- This 
parameter specifies the standard deviation of the wind 
azimuth angle in units of radians or degrees. If the 
value input is greater than or equal to "1", the program 
assumes the value is in degrees. If the value input is 


less than "1", the program assumes the value is in radians. 


Standard Deviation of the Wind Elevation Angle -- This 
parameter specifies the standard deviation of the wind 
elevation angle in units of radians or degrees. If the 
value input is greater than or equal to "1", the program 
assumes the value is in degrees. If the value input is 


less than "1", the program assumes the value is in radians. 


Time to Cloud Stabilization -- This parameter specifies 
the time to spray cloud stabilization in seconds. If 
this parameter is input as "0" or omitted from the input 


data, the program defaults TAU to 2.5 seconds. 


Measurement Time for SIGAP -- This parameter specifies 
the averaging time for the standard deviation of the wind 
azimuth angle in seconds. If this parameter is input as 
"O'" or omitted from the input data, the program defaults 


TAUO to 600 seconds. 
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Parameter 
Name 


SIGXYZ 


XLRZ 


DELU 


THETA 


DAREA 


Standard Deviation of the Spray Material Distribution -—- 
This parameter specifies the standard deviation of the 
spray material distribution along the spray line in units 
of meters. If this parameter is input as "0" or omitted 
from the input data, the program defaults SIGXYZ to 
WNGSPN/4.3. 


Lateral/Vertical Reference Distance -- This parameter 
specifies the lateral and vertical reference distance in 
meters downwind from the spray line where the spray cloud 
has dimensions given by SIGXYZ. This parameter is normally 
calculated by the FSCBG program. However, if XLRZ is 

input greater than or equal to zero, the input value is 


used. 


Wind Speed-Shear -- This parameter specifies the wind speed 
shear between the canopy and the aircraft height in meters 


per second. The default for this parameter is "0". 


Surface Mixing Layer Height -- This parameter specifies 


the height of the surface mixing layer in meters. 


Wind Direction -- This parameter specifies the wind 
direction (direction from which wind is blowing) in 


degrees, measured clockwise from 0 degrees (north). 


Area Assignment for Discrete Receptors -- This parameter 
specifies the area assignment in square meters to be used 
in the calculation of area-coverage for discrete receptors. 


If this parameter is set equal to "0" or omitted from 
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Parameter 


Name 
DAREA the input data, the program defaults DAREA to 10000 
(Cont. ) Square meters (1 hectare). 
BETAI1 Ratio of Lagrangian to Eulerian Time-Scales - This para- 


meter specifies the ratio of Lagrangian to Eulerian 
time-scales to be used in the calculation of the correction 
factor for SIGEP and SIGAP for crossing-trajectory effects 
of heavy drops when evaporation is not occurring. If 

this parameter is set equal to "0" or omitted from the 

input data, no correction factor is applied. Normally a 
value of "1" is input for the application of this correction 


factor. The equations used in the correction factor are 


given by 
stcep/|1+(BETAL-V,/ 
91°25 
(WSOCAN+SIGEP )) ; if BETAl > 0 & 
V,/ (SIGEP*WSOCAN) > .2386 
(non-evaporation only) 
SIGEPR = (3-11) 
SIGEP > df BETAL = 0c0r 


V,/(SIGEP*WSOCAN) < .2386 
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Parameter 
Name 


srcap/|1+4- (BETAL: V./ 


eee 
(WSOCAN* SIGEP)) ; if BETAl > 0 & 
(SIGEP*WSOCAN) > .2386 


(non-evaporation only) 
SIGAPR = (3-12) 


SIGAP ; of BETAl = 0 or 
VS/(SIGEP*WSOCAN) < .2386 


Receptor X (East-West) Coordinates -- This parameter is an 
array specifying the values of the X axis of the rectangular 
grid system and discrete receptor coordinates in meters. 

The first NXPNTS values input to this array comprise the X 
axis of the grid system and are input in ascending (negative 
to positive) order of value. These values are immediately 
followed by NXYPNT values giving the X coordinates of the 
discrete (arbitrarily placed) receptors. The negative X 
axis is west (270 degrees) and the positive X axis is 


east (90 degrees). 


Receptor Y (North-South) Coordinates -- This parameter is 

an array specifying the values of the Y axis of the rect- 
angular grid system and discrete receptor coordinates in 
meters. The first NYPNTS values input to this array comprise 
the Y axis of the grid system and are input in ascending 
(negative to positive) order of value. These values are 


immediately followed by NXYPNT values giving the Y coordinates 
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Parameter 
Name 


GAMA, GAMB 1 ; if V, > VSGAM(1) 
j = raat Sea 13 Pe > VSGAM(2) & Ve < VSGAM(1) 
(Cont. ) See lie ve < VSGAM(2) 


Three coefficients of GAMA, GAMB and GAMC each are required. 
The first value of each array is for settling velocities 
greater than or equal to the paramter VSGAM(1). The second 
value of each array is for settling velocities less than 
VSGAM(1) and greater than or equal to VSGAM(2). The third 
value of each array is for settling velocities less than 
VSGAM(2). Default values for these arrays are provided only 
if the parameter is omitted from the input data. The default 


values are 


GAMA(1) = 0.75 » GAMB(1) 
GAMA (2) 0.83465302, GAMB(2) 
GAMA (3) 0.91639996, GAMB(3) 


- 2.5 » GAMC(1) = 0.0 
- 6.9031391, GAMC(2) = 57.409256 
—22.357124 , GAMC(3) 821.42651 


VSGAM Settling Velocities for Fraction of Surface Reflection -- This 
array specifies two drop settling velocites in meters per 
second that break the curve of the fraction of material re- 
flected at the surface given by GAMA, GAMB and GAMC into three 
parts. Vv. values greater than or equal to VSGAM(1) use GAMA(1), 
GAM8(1) and GAMC(1) to calculate the fraction of material re- 
flected at the surface. Ve. values less than VSGAM(1) but 
greater than or equal to VSGAM(2) use GAMA(2), GAMB(2) and 
GAMC(2) to calculate the fraction of material reflected at the 
surface. ME values less than VSGAM(2) use GAMA(3), GAMB(3) 
and GAMC(3). Default values for VSGAM are provided only if 


VSGAM is omitted from the input data. The default values are 
VSGAM(1) = 0.04 , VSGAM(2) = 0.012 
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Parameter 
Name 


Y of the discrete (arbitrarily placed) receptors. The 
(Cont .) negative Y axis is south (180 degrees) and the positive 


Y axis is north (0 degrees). 


DX Source X (East-West) Coordinates -- This parameter is an 
array specifying the start and end X coordinate of each 
line source in meters. Input the start followed by the 
end X coordinate for each source. Either end of the line 


may be the start coordinate. 


DY Source Y (North-South) Coordinates -- This parameter is 
an array specifying the start and end Y coordinate of 
each line source in meters. Input the start followed by 


the end Y coordinate for each source. 


PCTMAT Fraction of Total Material -- This parameter is an array 
specifying the fraction of the total spray material in 
each drop-size category. Input these values in descending 


order of drop-size categories. 


GAMA, GAMB Coefficients of the Equation of the Fraction of Surface 
Reflection -- These three arrays specify the coefficients 
of the equation giving the fraction of material reflected 


at the surface as a function of the drop settling velocity 


R= CAMA, “© CAMB. ©.V-.4+'CaMC. - V7 (3-13) 
J aj s J Ss 


where 


V_= drop settling velocity in meters per second 
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Parameter 
Name 


DOSLV Dosage Area-Coverage -- This parameter is an array speci- 
fying the dosage levels in output dosage units for 
which dosage area-coverage is to be calculated. This array 
is not used by the program if either ISW(6) or ISW(12) are 


equal to '0". There are a maximum of 10 values possible. 


CONLV Concentration Area-Coverage -- This parameter is an array 
specifying the concentration levels in output concentration 
units for which concentration area-coverage is to be calcu- 
lated. This array is not used by the program if either 
ISW(7) or ISW(12) are equal to "0". There are a maximum of 


10 values possible. 


DEPLV Deposition Area-Coverage -- This parameter is an array speci- 
fying the deposition levels in output deposition units for 
which deposition area-coverage is to be calculated. This 
array is not used by the program if either ISW(8) or ISW(12) 


are equal to "0". There are a maximum of 10 values possible. 


Jokes Output Information 


The FSCBG program generates four categories of print output. 
Each category, except the first, is optional to the user. In the following 
paragraphs, each category of output is related to the specific input param- 
eter or parameters that control the output category. All program output 


is printed. 


Input Parameter Output. The first page of output produced by the 
FSCBG program is a listing of all input data read by the program. An example 
listing of this table is shown in Figure 3-1. Note that input parameters 
that have been omitted and are to be calculated by the program are shown as 


* symbols. A similar listing on a different computer might show $ symbols 
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rather than * symbols. Figure 3-1 consists of three pages of program out- 
put. The first page shows the parameters used by all program models, those 
used by the wake-settling velocity model and those used by the evaporation 
model. The second page of this table shows the parameters used by the 
canopy penetration model and the third page shows those used by the 


dispersion models. 


Evaporation Model Output. This output is controlled by the 
parameter ISW(5). The user has the option of not printing (ISW(5)="0" or 
printing the results of these calculations (ISW(5)="1"). An example out- 
put listing produced by ISW(5) equal to "1" is shown in Figure 3-2. The 
first page of Figure 3-2 shows calculated parameters used by the evapora- 
tion model that are primarily of theoretical interest. The second page 
shows a detail of the above canopy evaporation for drop size category l. 
The listing gives the height of the drop, the maximum, mean and minimum 
size, the total mass loss due to evaporation, the mass loss in the height 
interval due to evaporation, the travel time, the horizontal distance, the 
settling velocity and the percent of the material reaching the indicated 
height as drops. Following the evaporation detail for drop-size category 1, 
the program prints the coefficients for the equations of settling velocity 
as a function of distance, height as a function of distance, fraction 
of drops as a function of height, time as a function of height and drop 
size as a function of time. The evaporation detail and following equation 
coefficients are repeated for each drop-size category, largest to smallest. 
The program then prints the same information for each drop-size category 
below the canopy. The last three pages of Figure 3-2 show the evapora- 
tion detail below the canopy for drop size category 1. If the drop-size 
category evaporates prior to reaching the canopy the program will print 


a "DROP EVAPORATED" message for the below canopy evaporation detail. 


Canopy Penetration Model Output. This output is controlled by 
the parameter ISW(3). The user can elect either not to print (ISW(3)="0") 


or to print ISW(3)="1") the results of these calculations. An example output 


listing produced by ISW(3) equal to "1" is shown in Figure 3-3. The canopy 
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penetration model output gives the fraction of drop material that 
reaches each tenth of the forest canopy height and the horizontal dis- 
tance from the point of entry to the given tenth of canopy height. 

This information is printed for each drop size category. If a drop-size 
category evaporates prior to entering the canopy, the fraction of drop 


material in the table will be zero. 


Dispersion Model Output. This output is controlled by the 
parameter ISW(4). The user does not have the option of not printing these 
calculations. An example output listing produced by ISW(4) equal to "1" 
is shown in Figure 3-4. This figure shows a listing of the calculation 
of ground-level deposition. The listing gives the units, the maximum 
value and location of the maximum value followed by the deposition cal- 
culated at each receptor location. An output listing for dosage and/or 


concentration would have the same structure and content. 


Bo? USER'S INSTRUCTIONS 


367460 Program Description 


The FSCBG program is designed to calculate the wake-settling velo- 
city of the spray aircraft, the evaporation of spray drops, the percentage 
of spray material penetrating to fixed levels within a forest canopy and 
dosage, concentration and deposition at specified receptor coordinates 
downwind of multiple spray lines. The FSCBG computer program is written 
in FORTRAN IV and is designed for use on a UNIVAC 1108 computer. The 
program requires approximately 24000 words of executable core for instruc- 
tion and data storage on the UNIVAC 1108. The program requires data 
card image input (logical unit 5) and print output (logical unit 6). 


The program consists of a main program and 22 subroutines, (FSDTA, 
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CBGS1, CBGS2, CBGS3, EVPMD, EVPML, REGRS, AVGRD, IOPUT, UNRAN, CANPY, 
TEST, TESS, CBGS4, CBGS5, FOFX, ERFXS, VRICL, LNPN1, LNPN2, CBGS6 and 
TITLR). Although the program is written primarily for use on a UNIVAC 
1108 computer, the program is easily adaptable to other computers. The 
program can be segmented by using the subroutines beginning with CBGS 

and following non-CBGS subroutines as overlay segments. The program would 


then have 6 segments in addition to the main program segment. 


Be2 2 Control Language and Data Deck Setup 


Control Language Requirements. The following illustrates the 


required executive control statement runstream for a typical run on a 


UNIVAC 1108 operating system. 


Le @RUN,priority jobid,account ,userid,time,pages 


Optional, used to 
direct output to a 


Ze @SYM PRINTS, ,device ‘ i 
print device 


ee @ASG,A prog-file. 


Optional, used only 
when the FSCBG program 
input data has been 
placed in a file or 
data element within a 
file. 


4. @ASG,A data-file. 


Optional, used to di- 

rect print output to 

a specific device when 
running in demand mode 


on @ASG,CP print-file. 
@BRKPT  PRINTS/print-file 


6. @XQT prog-file.FSCBG 
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Us card-input-data 
or 
@ADD data 


or 
@ADD data-file.data-name 


8. @BRKPT PRINTS 
@FREE print-file. 
@SYM print-file,,device 


Je @FIN 


where 


priority 


jobid 


account 


userid 


time 


pages 


device 
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input FSCBG data on 
card or typed from 
a terminal 


FSCBG input data have 
been placed in a file 
or element within a 
file 


Optional, used with 
5 above to direct 
the print output to 
a specific print de- 
vice 


job run priority 


six-character user supplied job iden- 


tification 


account number 


up to 12 character user supplied pro- 


ject or user number 


execution time required in minutes 


output pages required 


printer symbiont name, onsite or re- 


mote, to which you desire the print 


file to go 


prog-file = the name of the program file. This 
illustration assumes the user has 
assembled and collected the FSCBG 
program into this file and called 


the absolute program FSCBG 


data-file the name of an optional data file 
into which the user has placed the 


input card image data for FSCBG 


i] 


print-file optional, user supplied, file name 
to be used for the FSCBG print output 


file. 


card-input-data FSCBG program input card image data 
defined in Section 3.1.2 and shown in 


Figure 3-5. 


Data Deck Setup. The card image input data required by the FSCBG 
program depends on the program options desired by the user. The card image 
input deck may be partitioned into five major groups of card data. Figure 3-5 
illustrates the input deck setup. Note that some of the card groups shown 
may be omitted from the input deck, depending on the program options chosen. 


The five major input groups are: 


ie Required Input data. These data consist of 1 or more data card 
images that are always included in the input deck. Figure 
3-5 shows the 5th and 6th card image of this group are in- 
put only if a specific parameter on the preceeding card is 


zero. 


Ds Wake Model Input Data. Only one of these two data card 


images is included in the input deck, depending on ISW(1). 
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included only if .VAPINF 
VRELHMO , RELHMU on the previous card is 
not input 


(1) 

included only if AIRDEN 
from the previous card 
of this group is not 
input or, "0" 


VAIRMOL, AIRPRS , VAPINF 


VIFWATR,WNGSPN,...,AIRDEN 


VISW( program control option array) 


VTITLE (array) 


VSQLST1 


(Required Input Data) 


FIGURE 3-5. Input data deck setup for the FSCBG program. (the V smybol 
means blank or no punch). 
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(2) 


VWAKVEL included only if ISW(1)="0" 


VARCRWT , ARCRSP included only if ISW(1)="1" 


(Wake Model Input Data) 


FIGURE 3-5. (Continued) 


(the Y symbol means blank or no punch) 
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VYDCL (array) 


VDBL (array) 


VDAL (array) include only if 


ISW(2)="2" 
VDCU (array) 


VDBU (array) 


YDAU (array) 


(3) 


include only if 
IsW(2)="1" and 
TFWATR="'2" 


VDFUSIV,HETLAT,CONMOL,...,CCONST ' 


VAIRPRS , AIRMOL ,VAPMOL, ... , RELHMU 


YDRPLWR (array) 


include only if 
ESW(2) <1"! 


include only if 


j Sr ees Isw(2)<"1" 
YDRPUPR (array) = 


(Evaporation Model Input Data) 


FIGURE 3-S. (Continued) 


(the V symbol means blank or no punch) 


Oi: 


VTREENV (array) 
(4) | YCOLEFF (array) 


YTREDEN (array) 
included only if 
: ISW(3) > Non 

VPRBPEN (array) 


VHGTCAN (array) 


VWSIN14 ,WSIN24 ,WSIN34 ,WSING4 


(Canopy Penetration Model Input Data) 


FIGURE 3-5. (Continued) 


(the VY symbol means blank or no punch) 


a2 


required end 
card 


Y SEND 


included only 
if, ESWiGS) > 0! 


VDEPLV (array) 


VCONLV (array) 


included only included only 
Lee eESW Ci Or icf SESWIGL2)i> 50" 


included only 


VYDOSLV (array) if ISW(6) >"0"" 


V GAMA,GAMB ,GAMC,VSGAM (array) 
VYPCTMAT (array) 


VDY (array) 


VDX (array) 


VY (array) 
included only 
if + TSWC4)2"0"" 
VX (array) 
V HM, THETA, DAREA, BETA1 


V SIGXYZ, XLRZ , DECAY , DELU 


VSIGAP, SIGEP , TAU, TAUO 


VNSOURC , NXPNTS ,NYPNTS , NXYPNT ,Q,SWATH,Z 


(Dispersion Model Input Data) 


FIGURE 3-5. (Continued) 
(the Y symbol means blank or no punch) 
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Or Evaporation Model Input Data. The data cards required in 
this section of the deck depend on the parameters ISW(2) 
and IFWATR. Two or more of these data cards will always 


be included in the input deck. 


4. Canopy Penetration Model Input Data. These data cards are 
included in the input deck only if option ISW(3) is greater 


than zero. 


De Dispersion Model Input Data. These data cards are included 
in the input deck only if option ISW(4) is greater than 
zero. Also, the last three cards of this group depend on 


ISW(12), ISW(6), ISW(7) and ISW(8). 
G. The last data card of the input deck must end with SEND. 
i eae Input Data Deck Format 


The format of the FSCBG program input data is the FORTRAN NAMELIST 
format and is described in detail in Section 3.2.4 below. The NAMELIST name 
is QLST1 and the first card of the input deck must have $QLST1 punched in 
columns 2 through 7 with columns 1 and 8 both blank. The FSCBG input param- 
eter names and values follow on the first and subsequent data cards and are 
terminated by SEND preceded by at least one blank. Column one of the input 
data deck must never be punched. Figure 3-5 shows the form of the input 
deck. These data may be entered directly into the runstream following the 
@XQT card or may be added (@ADD) from a previously created file or symbolic 
element within a previously created file. Table 3-1 gives the name of each 
input parameter, limits, default values and a summary description of each 
parameter. Also, Figure 3-6 gives a set of example input data coding forms 


for the FSCBG program. 
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TABLE 3-1 
FSCBG PROGRAM INPUT DATA PARAMETERS 


Limits Maximum 
Parameter Default * 


for Input Number of Description 


Rene Value Values Venus 


Optional Program Control Data 


PELE N/A oy Blanks Run Titling Information 
Characters 


ISsw(1) O51 il Wake-settling velocity 
option-—- 

O=use parameter WAKVEL 
for the wake settling 
velocity 

l=calculate the wake 
settling velocity 
from ARCRWT,WNGSPN, 
ARCRSP and AIRDEN 


TSW(2) Ora 2. Evaporation model option- 


O=spray drops are not 
evaporated, but used 
as input in all model 
calculations 

l=spray drops are evap- 
orated. Evaporated 
size is used in model 

calculations 

2=user inputs the coef- 
ficients of the equa- 
tions that describe 
drop size as a func- 
tion of time 


ISW(3) ORS 2 Canopy penetration model 
option- 


O=no canopy 

l=execute canopy pene- 
tration model 

2=execute and print 
calulations from can 


opy penetration 
model 


* if the parameter is omitted from the input deck, the parameter assumes 
its respective default value. 


oD) 


TABLE 3-1 (Continued) 


Maximum 


Parameter Heine eS 
for Input Number of 
Name 
Value 


Default * 


Description 


Dispersion models option- 


O=no dispersion model 
l=execute one or more 
of the dispersion 
models, depending on 

ISW(6),ISW(7) and 
ISW(8) 


Print evaporation model 
calculations option 


O=do not print calcu- 
lations 
l=print calculations 


Dosage model option- 


O=do not calculate 
dosage 
l=calculate dosage 


Concentration model op- 
tion- 
O=do not calculate con- 
centration 
l=calculate concentra- 
tion 


Deposition model option 
O=do not calculate 
deposition 
l=calculate deposition 


* If the parameter is omitted from the input deck, the parameter assumes 


its reactive default value. 
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TABLE 3-1 (Continued) 


Limits Maximum 
Parameter Default * feaee 
ans for Input Number of Description 


Value Values Values 


ISW(9) Dosage, concentration and 
deposition mass output 
units option- 

O=drops 
1=micrograms 
2=milligrams 
3=grams 
4=ounces 
5=pounds 


| ISW(10) Dosage time output units 
option- 
0=seconds 
1=minutes 


ISwW(11) Dosage, concentration and 
deposition volume or area 
length units output op- 
tion 

; area of 
3 
= ft 
= 3 acre 
= acre 
= hectare hectare 
ISsw(12) Onn Area-coverage option- 
O=do not calculate area 
coverage 
l=calculate area-cover- 
age for dosage, con- 
centration and/or de- 
position, depending 
on ISW(6),ISW(7) and/ 
or ISW(8) 


* If the parameter is omitted from the input deck, the parameter assumes 
its respective default value. 


D7. 


TABLE 3-1 (Continued) 


Limits ximu 
Ma 2 Default * 


for Input Number of Valtes Description 


Parameter 


ISW(13) 


ISW(14) 


ISW(15) 


ISW(16) 


* If the parameter is omitted from the input deck, the parameter assumes 
its respective default value. 
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TABLE 3-1 (Continued) 


Limits Maximum + 
Parameter for Input Number of VES Description 
Name Valse Values 
alu 


3 sig 


Isw(18) 


Isw(19) 


* If the parameter is omitted from the input deck, the parameter assumes 
its respective default value. 
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TABLE 3-1 (Continued) 


BAe Limits Maximum Deedee a, 
eee for Input Number of Description 


Value Values VET 


Required Input Parameters 


Spray liquid base 


l=water spray 
2=non-water spray 


Ee 


Aircraft wing span or 
helicopter rotor diam- 
eter (m) 


Wind speed above the 
canopy 


Zz 
~ 
> 


Height of aircraft above 
ground (m) 


Density of drop liquid 
(g em73) 


a 
rats 
ak 
a 


a1 
™. 
> 


Air temperature above the 
canopy (°C) 


AIRTPO Air temperature below the 
canopy (°C) 


% 


E 
A he 
ote ° 
> io) 


Air density (g ee 


Vv 
(2) 


Parameters that are Required if ISW(1)=0 


Wake-settling velocity 
(ms) 


Ei 


* If the parameter is omitted from the input deck, the parameter assumes 
its respective default value. 


** Jf this parameter is omitted from the input deck, the parameters AIRMOL, 


AIRPRS and VAPINF are required. Further, if VAPINF is omitted, the 
parameters RELHMO and RELHMU are required 
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TABLE 3-1 (Continued) 


Limits Maximum 
Default * nie 
for Input Number of Description 
Values 
Value Values 


Parameter 
Name 


Parameters that are Required if ISW(1)=1 


Aircraft ground-speed 


(ms ) 


Upper limit to the drop 
diameter for each drop 
size category. Input in 
descending order of size 
(micrometers) 


Lower limit to the drop 
diameter (micrometers) 
for each drop size cate- 
gory. Input in descend- 
ing order of drop size. 


ATRPRS 


Molecular weight of air 
ALRMOL (g mole~1) 


VAPMOL Molecular weight of evap- 
orating vapor (g mole~) 


Relative humidity above 
the canopy (%) 


Relative humidity below 
the canopy (%) 


* If the parameter is omitted from the input deck, the parameter assumes 
its respective default value. 


**K*k Default value is used if the parameter equals O or is omitted. 
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TABLE 3-1 (Continued) 


Limits Maximum 
for Input Number of 
Value Values 


Default * 
Values 


Parameter 


Description 
Name Pp 


DFUSIV Diffusivity of the evapo- 
rating vapor into air at 
the drop temperature 
(cm? s~) 


HETLAT Latent heat of vaporiza- 
tion (cal mole ~) 


Molal concentration of 
the air-liquid mixture 
(moles cnr 3) 


Thermal conductivity of 
the vapor into air at 
the drop Se eute 
—-{ oO 
cm K) 


(cal s 


Vapor pressure of the 
evaporating vapor at 
infinity (mb) 


BCONST, 21.07,*** | Vapor pressure equation 
CCONST 5249.9 constants 


* If the parameter is omitted from the input deck, the parameter assumes 
its respective default value. 


*xkk Default value is used if the parameter is O or omitted. 


xkkk The default value for these parmeters is calculated assuming similar to 
water. 


LOZ 


TABLE 3-1 (Continued) 


Maximum 
Number of 


Limits 
for Input Default * 


Values 


SUAS Description 


DAU, DBU, 
DCU 


Coefficients of the equa- 
tion of drop size as a 
function of time for each 
drop size category above 
the canopy 


D=DAU+DBU+ t+DCU> t* 


Input each array in de- 
scending order of drop 
size 


DAU ,DBU, 
DCU 


DAL, DBL, 
DCL 


Coefficients of the equa- 
tion of drop size as a 
function of time for each 
drop-size category below 
the canopy 

D=DAL+DBL* t+DCL« t“ 
Input each array in de- 
scending order of drop 
size 


WSIN44 


WSIN34 
se 


WSIN14 


WSIN34%*%% 


WSIN24%%% 


Wind speed in 
ter of canopy 


Wind speed in 
ter of canopy 


Wind speed in 
ter of canopy 


Wind speed in 
ter of canopy 


top quar- 
(ms ) 


third quar- 
(ms ~) 


second quar- 
(ms ~) 


bottom quar- 


(m s+) 


* if the parameter is omitted from the input deck, the parameter assumes 
its respective default value. 


*** Default value is used if the parameter is 0 or omitted. 
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Limits 
for Input 
Value 


co 


Parameter 


Name 


PRBPEN 


eal 


COLEFF >0 & <1 
or 
<0 

TREENV >0 


TABLE 3-1 (Continued) 


Maximum 
Number of 
Values 


me 


100 


Default * 
Values 


N/A 


ee i" 


N/A 


=13 


N/A 


Description 


Height of each forest 
canopy height class (m) 


Probability of drop pene- 
tration (fraction) for 
each canopy height class 


-1 
Tree density (stems acre ) 
for each canopy height 
class 


Collection efficiencies 
(fraction) for each drop 
size category 

or 
Diameters (cm) of the 
vegetative elements re- 
presentative of each 
canopy height class. 
Input as negative val- 
ues 


Tree envelope (width) (m) 
at each meter of tree 
height for each canopy 
height class from 1 meter 
above the ground to 
HGTCAN. Input HGTCAN(1) 
values, followed by 
HGTCAN(2) values, fol- 
lowed by HGTCAN(3) 

values 


Parameters that are Required if ISW(4)70 


NSOURC 


>0 & <100 


L&R 


Number of line sources or 
aircraft spray lines 


* If the parameter is omitted from the input deck, the parameter assumes 
its respective default value. 


**kk Default value is used if the parameter is O or omitted. 
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TABLE 3-1 (Continued) 


Peane cer Limits Maximum Deel # ee 
hee for Input Number of Description 


Value Values Values 


Number of X (east-west) 
coordinates in receptor 
grid system 


Number Y (north-south) 
coordinates in receptor 
grid system 


Number of discrete (arbi- 
trarily placed) receptors 


E 


Emission of spray material 
in gm if_$WATH-O, 

im gal vacre "if 

SWATH>0 


Distance between spray 
lines (m) if Q is in gal- 
lons/acre, otherwise omit 
or input as 0 


ert 
aake 


Height of dosage, concen- 
tration and deposition 
calculations (m) 


Standard deviation of the 
wind azimuth angle (radi- 
ans or degrees) 


pie 


Standard deviation of the 
wind elevation angle 
(radians or degrees) 


Time to spray cloud sta- 
bilization (s) 


Vv Vv Vv Iv 
(S) fo) ro) o> 
PS 


* If the parameter is omitted from the input deck, the parameter assumes 
its respective default value. 


*&*x*k Default value is used if the parameter is O or omitted. 


**** 100>NXPNISHNYPNTS+2*NXYPNT and 737>NXPNTS*NYPNTS+NXYPNT. 
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Parameter 
Name 


KKEEK 


TABLE 3-1 (Continued) 


Limits Maximum 
for Input Number of 
Value Values 


Default * 
Values 


KKK 
WNGSPN/4.3 


10000*** 


Tf the parameter is omitted from the imput deck, the parameter assumes 
its respective default value. 


Description 


Measurement time for 
SIGAP(s) 


‘Standard deviation of the 


source material along the 
spray line (m) 


Coefficient of time de- 
pendent exponential de- 
cay for the removal of 
material due to physical 
or chemical processes 


Lateral and vertical re- 
ference distance (m) 


Wind-speed shear_above 
the canopy (m s~1) 


Surface mixing layer 
height (m) 


Wind direction (degrees 
clockwise from north) 


Area assignment for dis- 
crete receptor points 
(m2) 


Default value if the parameter is O or omitted. 


This parameter is normally calculated by the program. However, if XLRZ 


is input > 0, the input value is used. 
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Parameter 


Name 


Limits 
for Input 
Value 


N/A 


N/A 


N/A 
N/A 


>0 & <1 


>0 


TABLE 3-1 (Continued) 


J Eesha Default * 
Number of 
Values 
Values 
a 
NXPNTS N/A 
+NXY PNT 


NYPNTS N/A 
4+NXYPNT 


at 
cate 
dell 


He 
o 


Description 


Ratio of Lagrangrian to 
Eulerian time-scales 
used in the correction 
factor on SIGEP and 
SIGAP for crossing-tra- 
jectory effects of heavy 
drops (fraction). Re- 
commended value is 1.0, 
if used 


X axis of receptor grid 
system followed by dis- 
crete receptor X coordi- 
nates. X axis should be 
in ascending order of 
value (m) 


Y axis of receptor grid 
system followed by dis- 
crete receptor coordi- 
mates. Y axis should be 
in ascending order of 
value (m) 


Start and end X coordi- 
nate of each source (m) 


Start and end Y coordi- 
nate of each source (m) 


Fraction of total spray 
material in each drop- 
size category 


Settling velocities that 
divide the curve of sur- 
face reflection coeffi- 
cient as a function of 
settling velocity into 
three curves (ms ~) 


* If the parameter is omitted from the input deck, the parameter assumes 
its respective default value. 


*** Default value if the parameter is O or omitted. 
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TABLE 3-1 (Continued) 


Limits 


Parameter 
Noe for Input 
Value 
GAMA , GAMB , N/A 
GAMC 


Maximum 
Number of 
Values 


each 


GAMA=.75, 


Default * 
Values 


-83465302, 
- 91639996 


GAMB=-2.5, 


=6. 9031391, 
=22 35/124 


GAMC=0.0, 


57.409256, 
821.42651 


Description 


Coefficients of equations 
using the surface reflec- 
tion coefficient as a 
function of drop settling 
velocity 


VS>VSGAM (1) uses GAMA(1), 
GAMB(1) & GAMC(1) 

VSGAM(1) >VS>VSGAM(2) 
uses GAMA(2), GAMB(2) 
& GAMC(2) 

VSGAM(2)>VS uses 
GAMA(3), GAMB(3) & 
GAMC (3) 


Parameters that are Required if ISW(12)>0 


DOSLV N/A 


10 


N/A 


: : : 
: ; 


* Tf the parameter is omitted from the input deck, the parameter assumes 


its respective default value. 
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N/A 


Dosage levels in units 
same as output dosage for 
| dosage area-coverage (re- 
quired only if ISW(6)>0) 


Concentration levels in 
units same as output con- 
centration for concentra- 
tion area-coverage (re- 
quired only if ISW(7)>0) 


Deposition levels in units 
same as output deposition 
for deposition area-cover- 
{age (required only if 
| ISW(8) 70) 


ee ee ek 
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B\o Ao NAMELIST Input Format 


The FSCBG computer program uses the FORTRAN NAMELIST DATA 
input format for input of control data. The namelist input data must be 
in a specific form in order to be read using a NAMELIST list. The first 
character (column 1) in each input card must be blank. The first card 
in the namelist deck contains the namelist name QLST1, preceded by the 
character $ (QLST1). The $ must be punched in column 2 and the name 
QLST1 must be followed by at least one blank. The last card in the 
input deck contains SEND, preceded by at least one blank, to terminate 
the control data input. Use the character & rather than $ on IBM computers. 


The form of the data items within the deck may be: 


a. Variable Name = Constant, -- The variable name may be a 
subscripted array name (ISW(4)) or a single variable name WAKVEL. Sub- 
scripts must be constants. The Constant may be an integer (whole 
number) or may be real (fractional part) and is terminated by a comma. 
Integer input parameter variables are those that begin with the letter 
I, J, K, L, M or N and must have a whole number Constant value. Real 
input parameters are those variables that begin with any other letter 
and may have a fractional part (45.342). Also, real parameters may have 
a Constant of the form X.XXEII, where E denotes the power of 10 of the 


value. 
Example: 


IFWATR=1,WNGSPN=9.144,WSOCAN=4.6, 
HGTCFT=30, 


b. Array Name = Set of Constants, -- The Array Name is not 


subscripted. The Set of Constants consists of constants of the type 
real or integer, depending on the Array Name, and all are separated by 


commas. The number of constants must be less than or equal to the 


maximum allowed by the program for the specific Array Name. Successive 


occurrence of the same constant value can be represented in the form 


k*¥Constant. 


Example: 


TSW=0,152,151,2*0,151, 
DRPUPR=414.0, 365, 341,320,304, 282,265,249, 


Cx Array Name = nnH ... Hollerith String ..., -- The Array Name 
is not subscripted (or is a simple Variable Name). The Hollerith String 
consists of labeling information or format information of any characters 
desired. There are exactly nn characters of information (including 
blanks) following the H and preceding a terminating comma. The only 


FSCBG input parameter that may use this form is TITLE. 


Example: 


TITLE=37HRUN TITLE FOR FSCBG PROGRAM EXECUTION, 


Input parameters and array elements that are not included 


(punched) in the input deck assume their respective default values. 


ore re Program Run Time and Page Output Estimates 


This section gives approximations to the computer run time and 
page output for the FSCBG program. Because of the variability of problem 
runs and input parameters, the equations in this section are meant only to 


give an approximation of the upper limit of the time and page usage function. 


a. Run Time. The run time required for a problem run is given 


by 
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Time (Seconds) 


—s 


Ny G o EL ae 38) ae Gey 6 Lae) 
Naa CGeve she Ne tof ))i47 (x See) 


yi 


oN ae (Cz. eH ,))| 


number of drop size categories 


humidity above the canopy (RELHMO) 


humidity below the canopy (RELHMU) 


where - 
Ng = the 
H = the 
fe) 
H = the 
u 
No = the number 
N. = the number 
N, = the number 
Be = the number 
N = the number 
xy 
N. = the number 
2D Boe? 
{1 = 
0 
1.6 
f., = 
0.4 


of 


of 


of 


of 


of 


of 


tree height classes (multi-storied) canopy 


trees (stems) per acre (TREDEN) 


receptors in the grid system x axis (NXPNTS) 


receptors in the grid system y axis (NYPNTS) 
discrete receptors (NXYPNT) 


sources (NSOURC) 


TE MUSWIC2Z i=l 


if TSWC2)Ae"1" 


LE -USWe2) =" 


Lf PuSWG2) 7-1 


Ly 


(3-14) 


aie 


Hh 
Ih 


ISW(2)="1" 


ISW(2)#"1" 


Isw(2)="1" 


IswW(3)="0" 


ISW(2)#"1" 


Isw(3)eV0! 


ISW(3)="0" 


TSW(3)>"'0" 


Isw(3)="0" 


or 


ISW(3)>"0" 


ISW(3)="0" 


TSwc3)e"0" 


Isw(3)>"0" 


The user should add any other time expected from other processing 


within the job to the estimate calculated via Equation 3-14. 


Also, if the 


values of f. above are not appropriate for your computer, recalculate and 


replace the f. values. 


De Page Output. 


FSCBG program depends on the problem being run and is given by 


Pages 


The number of pages of print output from the 


= A+B+CH+D 
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(3-15) 


where 


il Sn ecey SWS) = Ore 
IN = 
| 2et nN, 3 LE SWi(>)2°0": 
| 
| 1 SeieeiSWi(S)="OlvtorSWwiCs)=10" 
| B= {2+3-°N, 3; if ISw(5)>"0" & Isw(3)>"0" 
D ; if ISW(5)="0" & ISW(3)>"0" 
Y ; if Isw(4)="0" 
Cx = 
N 
5 + i -(9 +N ) 
it —=—___© | - (1sw(6)+Isw(7)+Isw(8)] ; if Isw(4)>"0" 
57 
0 ; If ISW(12)="0" 
D = 


ISW(6)+ISW(7)+ISW(8) ; if ISW(12)>"0" 


The user should add any other pages expected from the remainder 


of the job to the estimate calculated via Equation 3-15. 


* The[]symbols mean to round up to the next integer if there is any fractional 
part. 
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SECTION 4 
FSCBG EXAMPLE PROBLEM 


oak FSCBG EXAMPLE PROBLEM DESCRIPTION 


As we have indicated, the FSCBG model computer program can be 
used to calculate the deposition of drops within and below forest canopies 
for comparison with measurements. In this example problem, we use the 
dispersion and deposition models to calculate drop deposition along a 
line of deposition card samplers placed at ground level beneath a southern 
slash pine seed orchard canopy. A schematic diagram of the orchard and 
position of the card sampling line is shown in Figure 4-1. The rows of 
trees are separated by 30 feet and the spacing between stems along each 
row is 15 feet. The spray aircraft, a Stearman biplane, flies north and 
south between rows of trees beginning between rows 11 and 12 and ending 
between rows 6 and 7. The aircraft is assumed to fly at an altitude of 
60 feet above ground and to spray a water-based pesticide at a rate of 
57.33 grams per meter. The rectangular coordinate system used in defining 
the problem has its origin at the first tree in row 1 shown at the lower 


right-hand corner of Figure 4-1. 
Meteorological Input Parameters 


For the example problem, which is based in part on an actual 
experiment conducted by the Forest Service, we have assumed that meteoro- 
logical measurements were made above the canopy using a bidirectional wind 
sensor and a temperature sensor mounted on a 18-meter tower. An analysis 
of the recorded wind azimuth and elevation data over a 10-minute measure- 
ment period (T=600 seconds) beginning at the time the aircraft started to 
release the spray material yielded estimates of the standard deviations 
of the wind azimuth angle (O,) and wind elevation angle (9,) shown in 


Table 4-1. The mean wind direction (9) and wind speed (u) shown in 
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Example spray problem design. 


FIGURE 4-1. 
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TABLE 4-1 
METEOROLOGICAL INPUT PARAMETERS 


0, {tT =600s} (deg) SIGAP 
On (deg) SIGEP 
Ho (m) HM 
= = 
ie ame: Se) WS@CAN 
8 (deg) THETA | 
- aT 
Nu Guss: =) DELU 
- (mb) AIRPRS 
RE 3CZ) RELHMG 
Oo 
i Co) AIRTP@ 
u (m Bas) WSINI4 
esl 
u Ger) WSIN24 
4632 as 
=i 
1.33 (ns) WSIN34 
(m s+) WSING4 
Yos4 co 
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Table 4-1 were also estimated for the 10-minute period from the wind 
sensor. The wind speed shear (Au) above the canopy was set equal to zero 
because wind speed was measured only at a single height and furthermore 
the shear was not expected to have a significant influence on the calcu- 
lations. The temperature above the canopy (T,) was 14.4 degrees Celsius. 
The relative humidity (RH) was measured near ground level. The measured 
values of air temperature and relative humidity shown in Table 4-1 were 
assumed in the example problem to represent the temperature and relative 
humidity both above and within the canopy. The wind speeds within the 
canopy (a, shown in the table are based on a normalized height profile 


for winds in an isolated conifer stand presented by Fritschen, et al. 


COTO). 


Source Input Parameters 


Source input parameters used in the model calculations are shown 
in Table 4-2. The source dimension O, was determined form observations that 
the spray cloud quickly attained a dimension equal to about 1.5 wingspans 
behind the aircraft. The value of O, was determined by dividing this 
dimension by 4.3 (9 =1.5b/4.3) under the assumptions that the distribution 
of material within the spray cloud was Gaussian and that the visible edge 
of the cloud represented the point at which the concentration was one- 
tenth the concentration at the cloud centroid. The aircraft air speed 
was approximately 90 miles per hour (40.2 m SED along the flight lines 
shown in Figure 4-1. The flight lines were longer than shown in the 
schematic diagram, each line actually extending from 384.3 meters south 


(y=-3847.3 m)! t0,3/70 meters north (y=370.3) of the first tree in Row 1. 


Table 4-3 shows the spray drop-size distribution used in the 
example calculations. This distribution is based on drop-size measurements 
made during trials in which the Stearman aircraft disseminated material 


over a sampling line in an open area. 
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TABLE 4-2 
SOURCE INPUT PARAMETERS 


HGTCFT 
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Drop Upper 
Limit ( m) 
DRPUPR 
1420 
1020 
840 
742 
667 
Di 
514 
456 
3a 
302 
253 
200 
169 
135 
100 


51526 


TABLE 4-3 
SPRAY DISTRIBUTION PARAMETERS 


Drop Lower 
Limite ‘© m) 
DRPLWR 
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Fraction of Material 
In Drop-Size Category 
PCTMAT 
0.01 
0.02 
(B}A0)8) 

0.04 
O86 
0.10 
0.10 
0.20 
0.10 
Ol 
0.10 
0.04 
0.03 
0.02 
0.009 


0.001 


Forest Description Parameters 


Forest description parameters for the example problem are shown 
in Table 4-4. Since the spacing between trees is 15 feet along the rows 


and there is 30 feet between rows, the tree density is 


1 stem - 43560 £t2 a 96.8 stems 


(30x15) Pen acre acre 


The average measured height of the trees (HD it the: vicinity of the 
sampling line was 12 meters. The tree widths W, for heights from 1 to 
12 meters were obtained by averaging measurements made of four typical 
trees near the sampling line. In the calculations, the program was set 
to calculate impaction efficiencies ce) from Equation (2-77) using a 


default value of s equal to 13 centimeters. 


Ob ee FSCBG EXAMPLE PROBLEM INPUT DATA 


The input data used for this example problem are shown in Figure 
4-2 coded on example coding forms identical to those shown in Section 3. 
Parameters and arrays on the coding forms that are crossed out indicate 
parameters that are not required (WAKVEL), were previously specified 
(RELHMO on sheet 2), are being defaulted (AIRTPU) or are trailing 
comma's that have no associated coded value. The first data card image 
begins with SQLST1 (the name of the NAMELIST) and the deck ends with 
SEND. The second input card shows the 47 character title. Following 
the title the program option switch ISW is coded. The option switch 
values shown in Figure 4-2 specify that the program is to calculate the 
wake-settling velocity, that spray drops are to be evaporated, a canopy 
is used and print output from the canopy model is desired, dispersion 


modeling is to be done, evaporation model results are to be printed, 


126 


TABLE 4-4 
FOREST DESCRIPTION PARAMETERS 


TREDEN 


D. (stems acre -) 


HGTCAN 


i=1 TREENV Oe 
2 O29 
3 Osis. 
4 3-8 
5 6.6 
6 Gi2 
7 3-8 
8 6.4 
9 5.0 
4.0 
303 
156 
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dosage is not calculated, concentration is not calculated, deposition is 
to be calculated, output mass units are drops, dosage time units do not 
apply, length output units are meters and, finally, area coverage is not 
calculated. Following the program options are the remainder of the 
"required" input parameters. A "required" parameter means that a value 
for the parameter must be either input or defaulted and does not mean 
that punching a value for the parameter is required. For example, the 
parameters DENLIQ, AIRTPU and AIRDEN have been omitted (crossed out) 
because a default value will be used for these parameters. Input parameters 
for the evaporation model and canopy penetration model are given on 
sheets (upper right corner) 1 through 3 of Figure 4-2. Those for the 
dispersion models are given on pages 4 through 5 and on the extra sheet 
following sheet 5. This extra sheet is provided for those arrays that 
exceed the space provided on sheets 1 through 5. For example, there is 
not enough room to code all of the values for the x axis of the receptor 
grid system on sheet 4, so the remainder of the x axis values have been 
coded on the extra sheet shown after sheet 5. Any values coded on this 
extra sheet that are continuations of values previously shown must be 
merged back into the input deck in their proper position prior to 


execution of the program. 


Ao FSCBG EXAMPLE PROBLEM OUTPUT 


The FSCBG example problem output is shown in Figure 4-3. Only 
selected pages of the program output have been included in Figure 4-3 
due to the size of the output listing (114 total pages). The first page 
of Figure 4-3 gives the FSCBG required input data, wake-settling velocity 
model input data and the evaporation model input data. Note that some 
parameters have *****k symbols for a value. This means the value was not 
specified on input and is to be defaulted or calculated by the FSCBG 


program. The second and third pages of Figure 4-3 show the input data 
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used by the canopy penetration model and dispersion models, respectively. 
The ***** symbols are also used on these pages to show parameters that 
are to be calculated by the program. Pages 4 through 111 are printed 
only if the parameter ISW(5) is set non-zero. In this example, only 
pages 4 through 7 and 63 through 66 (upper right corner) are shown. 

Page 4 gives calculated parameters used by the evaporation model. Pages 
5 through 7 show the above canopy evaporation calculations for drop-size 
category 1. A set of pages similar to pages 5 through 7 are printed for 
each drop size category for above canopy evaporation. The program 
prints the height, the maximum, mean and minimum drop diameter, the 
total mass lost to evaporation, the mass lost to evaporation in the 
height interval, the fall time, the horizontal travel distance, the 
settling velocity and the percent of spray material reaching the given 
height as drops. This information is printed down the page until the 
minimum drop diameter of 2 micrometers is achieved, or the drop is 
essentially evaporated or until the height falls to -10000 meters. The 
evaporation is calculated until these criteria are met because the FSCBG 
program requires the complete curve of evaporation in order to calculate 
the equations of height versus distance, drop size versus distance, etc. 
above the canopy. The last page for each drop size category (page 7) 
details the equations required by the canopy penetration model and 
dispersion models. Pages 63 through 66 (upper right corner) are the 
same in content as pages 5 through 7, except these pages give the evap- 
oration detail below the canopy for drop-size category 1. The below- 
canopy evaporation calculations are stopped at the ground and are not 
made for negative heights as in the case for the above-canopy evaperation 
calculations. Page 112 (upper right corner) of Figure 4-3 shows the 
output from the canopy penetration model and is obtained by the ISW(3) 
option. For each drop-size category, the FSCBG program prints the 
fraction of drops reaching each 10th of the forest canopy height. The 
fraction of drop material printed here accounts for losses due to evap- 
oration and losses due to impaction in the upper areas of the canopy. 


Note that drop size category 16 only reaches the 10.8 meter level. This 


iESyb 


is because the drop evaporates before it can penetrate any deeper into 
the canopy. Also, the horizontal travel distance required for drops in 
each category to reach each 10th of the forest canopy height is printed. 
Pages 113 and 114 (upper right corner) of Figure 4-3 show the ground 

level deposition calculated for this example problem. The output dis- 
plays "DEPOSITION", the output units, the maximum calculated value and 
location followed by the deposition at each input receptor point. A 


listing of dosage or concentration would appear in the same format. 


The FSCBG program has also been used to calculate deposition 
occurring along the sampling line shown in Figure 4-1 for relative 
humidities of 70 and 50 percent, with all other model inputs described 
in Section 4.1 remaining the same. The results of the ground-level 


deposition for all three humidity values are shown in Figure 4-4. 
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SECTION 5 
FSCBG PROGRAM LISTING 


This section gives a complete FORTRAN listing of the FSCBG com-— 


puter program and all subroutines, 
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DISCLAIMER STATEMENTS 


The use of trade, firm, or corporation names is for the information 
and convenience of the reader. Such use does not constitute an official 
evaluation, conclusion, recommendation, endorsement, or approval of any 
product or service to the exclusion of others which may be suitable. 


This publication reports research involving pesticides. It does not 
contain recommendations for their use, nor does it imply that the uses 
discussed here have been registered. All uses of pesticides must be 
registered by appropriate state and/or federal agencies before they can 
be recommended. 


